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ABSTRACT 
An  airborne  direction- finding  technique  capable  of 
handling  multiple  emitters  was  developed.   An  adaptive  gate 
size  was  introduced  in  the  track  correlations  based  on  co- 
variance  relations  of  established  tracks  and  observations. 
Track  files  were  generated  based  upon  the  existence  of 
various  target  parameters;  i.e.,  frequency,  pulse  repetition 
frequency,,  pulse  width  and  direction-of-arrival .   To  test 
the  angular  resolution  capability  of  the  filter,  emitters 
in  close  proximity  to  each  other  with  identical  electronic 
characteristics  were  used  in  the  simulation.   Target  loca- 
tions are  calculated  in  a  cartesian   coordinate  system 
where  the  sphericity  of  the  earth  is  taken  into  account, 
and  with  appropriate  coordinate  transformation  computational 
simplicity  is  preserved. 
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I.   INTRODUCTION 

In  airborne  direction-finding  (DF)  a  good  procedure 
for  locating  the  emitter  given  angle  observations  only  is 
of  primary  importance.   This  becomes  very  difficult  in  a 
multi-emitter  environment.   With  sophisticated  hardware  and 
software,  many  of  these  emitters  can  be  identified  and  sorted 
by  their  pulse  repetition  frequency  (PRF) ,  and  pulse  width 
(PW) .   However,  without  accurate  direction  of  arrival  in- 
formation and  methods  for  correlating  these  with  successive 
past  emissions,  the  problem  is  not  solved.   Central  to  all 
of  this  is  the  estimate  of  the  location  of  the  emitters. 

In  the  emitter  locating  portion  of  the  study,  it  was 
assumed  that  all  emitters  have  the  same  frequency,  PRF,  and 
pulse  width,  and  thus,  except  for  their  angle  of  arrival, 
all  correlate  as  the  same  target.   This  is  really  a  very 
difficult  problem.   To  get  an  idea  just  how  difficult, 
consider  a  rather  typical  environment  of  six  emitters  in 
reasonable  proximity,  and  an  aircraft  with  a  noisy  antenna 
system  making  ten  observations  of  each  of  these  emitters. 
The  possible  number  of  locations  with  a  straight-forward 
method  of  locating  by  calculating  all  possible  intersections 
of  bearing  lines  will  give  1,620  locations  computed  by 
the  formula 

2 

m  n (n-l)/2 

where  m  =  number  of  emitters  and  n  =  number  of  observation 
points.   The  computer  does  not  know  whether  there  are  six 


emitters  or  1,620.   Without  sophisticated  association  and 
correlation  routines  for  correlating  new  observations  with 
previous  tracks  and  adaptive  angle  tracking  and  gating, 
the  problem  is  virtually  insolvable. 

Two  methods  of  solving  the  location  and  correlation 
problem  were  investigated.   The  first  uses  a  Kalman  filter 
to  compute  predicted  values  of  angle  of  arrival  based  on 
past  data,  and  a  correlation  routine  to  correlate  the 
present  observations  with  these  predicted  values.   If  no 
correlation  was  possible,  it  initializes  a  new  track.   A 
correlated  observation  is  then  used  to  compute  an  optimal 
estimate  of  the  present  angle-of-arrival  for  that  track. 
In  the  second  method,  all  data  from  a  number  of  observation 
stations  are  processed  simultaneously  by  first  calculating 
all  possible  target  locations  based  on  intersections  with 
the  first  bearing  obtained.   Then,  the  mean  position  of 
this  cluster  is  found,  declared  a  probable  target  location, 
and  the  observations  having  contributed  to  this  location 
are  correlated  to  this  target.   In  both  cases,  the  target 
location  is  found  in  an  earth-centered  coordinate  system 
(where  the  x_  axis  is  the  polar  axis  and  the  x, ,  x~  plane 
is  the  equatorial  plane  of  the  earth,  with  the  positive  x, 
axis  corresponding  to  0°  longitude  and  the  positive  x~ 
axis  corresponding  to  90°  longitude  east) .   For  each  obser- 
vation, an  observation  plane  is  defined,  passing  through 
the  earth's  center  and  the  bearing  line  at  the  observation 
point,  such  that  in  the  noiseless  case,  this  plane  would 


also  pass  through  the  target  point.   The  i  ,  plane  is  then 
defined  by 

aliXl  +  a2iX2  +  a3iX3  =  °  '  (1) 

where  the  coefficients  a.,  can  be  found  from  the  observation 
station  coordinates  and  the  observation  angle.   The  procedure 
for  obtaining  a  target  location  is  to  take  the  equations  of 
the  observation  planes,  two  at  a  time,  with  the  constraint 
equation  that  the  emitter  be  located  at  the  earth's  surface. 
Solving  the  equations  of  two  planes  below,  gives  the  desired 
target  location. 

allxl  +  a12x2  =  "a13C  (2' 

a21Xl  +  a22x2  =  "a23°  (3) 

The  equations  were  modified  to  solve  for  x,  and  x-  with 
X-  —  c.   The  constraint  equation  is  accounted  for  by 
solving  for 

E  =  R/(x'2  +  x'2  +  c2)1/2  (4) 

where  R  =  radius  of  the  earth.   The  resulting  solution 

for  x.  is 

1 

x.  =  E  x!  i  =  1,  2,  3.  (5) 

By  setting  x~  =  c,  the  observation  planes  have  been  projected 
into  lines  on  the  plane  x,  =  c,  then  in  solving  (2)  and  (3) 
the  intersection  of  these  lines  in  the  x_  =  c  plane  is 
found.   The  parameter  E  is  then  utilized  to  project  this 
point  to  the  surface  of  the  earth.   With  proper  coordinate 
rotation,  it  is  possible  to  have  the  projection  plane  (x^=c) 
earth  surface  tangent  and  aircraft  centered. 
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II.       COMPUTATIONAL    PROCEDURE 

A.  VARIOUS  APPROACHES 

As  is  always  the  case ,  there  are  many  ways  to  attack 
this  problem.   Possibly  the  best  technique  for  evaluating 
the  various  alternatives  is  through  a  realistic  statistical 
evaluation  of  a  multi-emitter  environment.   But  first  a 
careful  description  of  the  two  alternative  techniques 
considered  here,  and  the  methods  of  analyzing  the  tech- 
niques (i.e.,  the  simulation  scenario)  is  in  order.   First, 
a  procedure  for  generating  the  target  and  receiver  kinema- 
tics is  necessary.   Then,  for  the  multi-emitter  environment, 
a  very  sophisticated  association  and  correlation  routine 
for  relating  the  observations  with  existing  tracks  is 
required.   Lastly,  the  location  of  the  emitters  must  be 
found,  given  either  estimates  of  bearing  angles  or,  alterna- 
tively, a  large  set  of  noisy  bearing  lines  giving  a  set  of 
noisy  locations  from  which  a  location  estimate  is  derived. 

B.  AN  ADAPTIVE  BEARING  ANGLE  FILTER 

A  Kalman  filter [1]  is  used  to  obtain  estimates  and 
prediction  of  bearing  angles  for  the  next  observation  point 
to  allow  the  correlation  of  observations  with  past  data. 
Each  observation  and  bearing  line  will  have  an  unknown 
angular  rate  associated  to  it.   This  rate  depends  on  the 
speed  and  heading  of  the  aircraft,  the  bearing  itself, 


and  the  unknown  range  to  the  emitter  for  fixed  or  slowly 
moving  targets.   It  is,  however,  approximately  known; 
and  this  approximation  is  used  in  the  initialization  of 
the  filter  for  each  target.   The  initial  uncertainity  of 
angle  and  angular  rate  on  track  initialization  is  accounted 
for  through  the  apriori  covariance  of  error  initialization 
in  the  filter.   The  filter  also  takes  into  account  missing 
data  points.   That  is,  at  any  instant  of  look,  at  most 
only  50%  of  the  targets  might  be  emitting.   Trackfile 
management,  like  initializing  a  new  track,  dropping  a  track, 
or  flagging  a  missing  observation  are  done  in  the  correla- 
tion section. 

The  noisy  observation   z   is  assumed  to  be  described 
by 


zk  =  hk_^  +  vk  ,       k  =  1,  2,  ...  (6) 


m 

Here  h,  =  (1   0),  and  Q^  =    (ekek)   is  the  true  observation 
angle  and  angular  rate,  and  vk  represents  a  white  noise 
sequence.   Since  the  observer  is  moving  in  time  (the 
subscript  k   denotes  the   k    instant  in  time) ,  a  dif- 
ference equation  is  used  to  approximate  the  system  dynamics 
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%  -  •k,k-l§k-l  * 


k  -  1,2,...  (7) 


No  random  forcing  is  included;  this  implies  stationary 
targets  and  exact  knowledge  of  the  observer's  position. 
The  filter  can  be  extended  to  take  moving  targets  into 
account  by  augmenting  the  difference  equation  above,  and 
using  an  appropriate  non-zero  Q  matrix  when  initializing 
the  recursive  filter  gain  and  covariance  equations.   The 
mean  and  variance  of  the  measurement  noise  sequence  v,  is 
given  by: 

E(vk)  =  0  for  all  k  (8) 

and       E^vkvj)  =  Rk6jk  *  (9) 

where      « ..  =  °'  j^ 
Dk    1,  j-k 

The  recursion  equations  for  filtering  are  defined  below, 

where  8.,.  denotes  the  estimate  of  the  state  ©_.  based  upon 

measurement  data  (z, ,  z», . . . . .  z.) 

1    2        j 

-k/k-1  =  $k,k-A-l/k-l  (10) 


Vk  =  i/k-1  "  2k(VUk0kA-l)  '  (11) 

where  g,  are  the  filter  gains. 

The  recursive  equations  for  the  filter  gain  and 
covariance  matrixes  are  defined  bv: 
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2k  "  pkA-l  £   (!ik  pkA-l  fik  +  V"1  (12) 


pk/k  =  Pk/k-l  "  %&pkA-l  (13) 

Pk+1A  =  *k+l,k  pkA  *k+i,k  +  Q  (14) 

Q  is  a  measure  of  random  exitation  in  the  process  and 
for  a  stationary  emitter  is  omitted. 

No  estimate  of  the  apriori  value  of  Q_   is  available. 
Therefore,  the  filter  is  initialized  with  the  first  measure- 
ment in  the  following  fashion: 

ex  =  zx  (15) 

0j_  =  (v/R)  sin(z1)  .  (16) 

In  equation  (16),  v  is  the  velocity  of  the  observer, 
R      the  distance  to  the  target,  and  z,  the  measured  bearing 
angle.   The  only  unknown  in  this  equation  is  R,  but  this  can 
be  approximated.   A  value  of  150  nmi  is  assumed  in  the 
computer  simulation. 

During  simulation  it  was  found  that  a  simple  $,  v-i 

1    f 


matrix  of  the  form   $,  ,  , 


does  not  adequately 


0    1 

represent  the  system.   This  can  be  seen  in  Figs.  2  and  3 
where  the  mean  errors  of  bearing  angle  and  angular  rate 
are  shown  using  the  simple  and  an  augmented  $  .  .  _-.  matrix 
in  a  Monte  Carlo  simulation.   The  augmentation  was  found 
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R  =I50NM 


=  360knts 


Fig,  1,  Simulation  scenario  at  first  observation  point  with 
three  of  six  targets  randomly  emitting  at  an  average  range 

OF  150  NMI . 

Initialization  equation:   9  =  (v/R)  sirr(9) 
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Fig.  2.  Mean  error  of  bearing  angle  obtained  from  a 
Monte  Carlo  simulation  on  the  Kalman  filter  with  a  single 
target  and  perfect  correlation  for  a  total  OF  20  STATION  POINTS. 

d  mean.  error  of  the  measurements 

+  mean  error  of  angle  estimates  using  a  nonlinear  ±  matrix 
in  the  Kalman  filter  prediction  equation 

x  mean  error  of  angle  estimate  using  a  linear  ±  matrix  in 
the  Kalman  filter  prediction  equation 
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Fig.  3,  Mean  estimate  of  angular  rate  and  mean  rate  error 


a   true  angular  rate 

+   mean  estimate  of  angular  rate  and  mean  rate  error 
obtained  when  using  nonlinear  ±  matrix 

x   mean  estimate  of  angular  rate  and  mean  rate  error 
obtained  when  using  linear  ±  matrix 
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by  introducing  an  additional  difference  equation: 

°k/k-l  =  VlA-l  +  ekA-lT  (17) 

The  addition  of  another  variable,  the  second  derivative  of 
the  bearing  angle,  does  not  greatly  complicate  the  computa- 
tion;  because  for  stationary  targets  (see  Fig.  1)  ,  8  can  be 
expressed  in  terms  of  already  defined  variables.   The 
augmented  recursion  equation  for  0ic/lc_-i  then  becomes: 


Vk-i =  ViA-i  (1+2*\-i/k-i*cot(9k/k-i)T)  <18) 


This  did  improve  the  filtering  results  and  eliminated  the 
bias . 

The  easiest  and  fastest  way  to  compute  an  apriori 
target  location  using  the  results  obtained  by  the  Kalman 
filter  is  to  compute  the  target  location  based  on  the  first- 
observation  made  and  on  the  estimate  of  bearing  angle 
obtained  at  a  later  observation  site,  after  a  number  of 
observations  have  been  made.   However,  the  first  observa- 
tion is  an  unfiltered  measurement,  and  therefore  may  have 
quite  a  large  error  which  will  degrade  the  apriori  target 
position  seriously.   To  overcome  this  difficulty,  the 
measurements  are  not  only  used  in  the  Kalman  filter  to 
compute  a  current  estimate,  but  also  to  obtain  a  smoothed 
estimate  of  the  first  observation.   This  smoothing  is  done 
by  a  technique  proposed  by  H.  E.  Rauch[2].   The  smoothed 
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estimate  of  the  initial  observation  after  k  measurements 
are  taken  is  given  by: 


h/k  =   °1A-1  +  Dl/k  %  (zk  "hk  9k/k-l) 


(19) 


T       -1 
Dl/k   Dl/k-l  Pk-l/k-l  $k,k-l  Pk/k-l 


(20) 


Equation  (20)  shows  that  if  the  simple  $  matrix  is  used, 
D,  «  is  data  independent  and  may  be  precomputed  and  stored 
For  this  case,  and  using  the  initial  values  given  in 
Appendix  A, 


'1A 


1    -(k-1) 
0       1 


and  since  ony  a  smoothed  value  of  the  bearing  angle  and  not 
the  angular  rate  is  required,  only  a  single  scalar  equation 
is  needed  to  obtain  a  smoothed  estimate  of  the  initial 
bearing  angle. 

This  leads  to  the  following  set  of  recursive  filtering 
and  smoothing  equations  in  scalar  form: 


6k/k-l  "  6k-l/k-l  +  Tek-l/k-l 


ek/k-l  ~9k-l/k-l  (1+2T8k-l/k-l  cot  ^k/k-1^ 


»k/k   =  6k/k-l  +  gk(1)  (zk  "  ek/k-l) 


(21) 


'k/k 


'k/k 


-1  +  gk(2)  (zk  "  ek/k-l] 


'i/k   "  Vk-i  +  5k(D(zk-ek/k_i)  -  (fc-^^^rVk-i* 
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C.   PARAMETER  CORRELATION 

The  purpose  of  correlation  is  to  assign  observations 
to  already  established  tracks  if  possible,  or  flag  an  obser- 
vation as  a  new  target  if  this  observation  does  not  correlate 
to  previously  established  tracks.   To  each  observation  (or 
track)  five  parameters  are  associated: 

Frequency 

PRF  (Pulse  Repetition  Frequency) 

PW   (Pulse  Width) 

Correlation  counter 

DOA  (Direction  of  Arrival) 
For  each  of  these  paramters  a  gate  is  selected.   The  size 
of  the  gate  depends  on  the  expected  accuracy  of  the  measure- 
ment and  the  confidence  associated  with  the  track  parameter 
in  question. 

First,  an  observation  is  correlated  to  a  track  if  the 
values  of  all  five  parameters  are  within  their  respective 
gate  width.   This  may  correlate  one  observation  with  several 
tracks,  or  several  observations  with  the  same  track.   To 
eliminate  this  multiple  correlation,  four  rejection  rules 
and  one  control  rule  are  used: 

1.  A  track  which  correlates  with  several  observations 
rejects  any  observation  held  in  common  with  another 
track  if  the  common  observation  is  the  only  obser- 
vation correlating  with  the  other  track. 

2.  A  track  correlating  with  several  observations, 

some  of  which  are  not  held  in  common  with  other 

tracks,  rejects  observations  held  in  common  with 

other  tracks. 
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3.  When  several  observations  correlate  with  one  track, 
the  closest  observation  is  correlated  to  that  track, 

4.  VThen  several  tracks  correlate  with  one  observation, 
the  observation  is  correlated  with  the  closest 
track . 

5.  If  an  observation  is  found  uncorrelated,  and  at 
the  same  time  there  are  cne  or  more  tracks  uncorre- 
lated  to  any  observation,  then  this  observation 

is  correlated  to  the  first  uncorrelated  track 

found  within  the  gates. 
The  gates  are  chosen  to  be  two  sigma  of  the  target  parame- 
ter measurement  noise  plus  two  sigma  of  the  confidence  inter- 
val associated  with  the  track  parameter  estimate  error 
which  is  a  function  of  the  number  of  observations  already 
correlated  to  that  track. 

Rule  5  was  included,  because  rules  3  and  4  can  reject 
too  many  of  the  initial  correlations.   If,  for  instance, 
two  observations  are  correlated  to  the  same  two  tracks, 
it  is  possible  that  one  of  these  observations--say ,  the 
first — is  closer  to  both  tracks  than  the  second  observation. 
Then  rule  3  rejects  the  second  observation  from  both  tracks 
and  leaves  the  first  correlated  to  both  tracks.   Rule  4 
then  selects  the  closest  track  as  correlated  to  the 
first  observation.   This  leaves  unconnected  the  second 
observation  and  one  track  which  was  correlated  initially. 


19 


D.   TARGET  LOCATION  AND  CORRELATION  BY  CLUSTER  POINTS 

The  intent  of  this  location  method  is  to  find  the  most 
probable  target  locations  based  on  observations  made  at 
several  stations,  flag  the  observations  leading  to  this 
location,  and  give  a  standard  deviation  to  reflect  the  quality 
of  the  obtained  observation.   Fig.  4a  shows  four  observation 
stations  (Points  I,  II,  III,  IV)  and  possible  bearings  at 
these  stations.   A  look  at  this  Figure  indicates  quite 
clearly  that  points  X,  Y,  and  Z  are  not  only  the  most 
probable,  but  also  the  only  logical  choices  for  target 
locations.   To  find  these  locations  by  means  of  a  digital 
computer,  one  observation  from  station  I  is  selected 
(Fig.  4b) ,  and  crossing  points  with  all  bearings  from  all 
other  stations  are  calculated.   Then,  the  mean  and  standard 
deviation  of  these  locations  is  found.   Since  in  the  close 
neighborhood  of  point  X,  three  points  are  located;  the  mean 
position  will  be  relatively  close  to  this  point,  but  the 
far  away  points  will  make  the  standard  deviation  large. 
Now  all  points  outside  a  1.1a  interval  from  the  mean  are 
eliminated,  and  with  the  remaining  points  a  new  mean  and 
standard  deviation  is  calculated.   This  procedure  will  elimi- 
nate erroneous  bearing  lines  and  will  finally  produce  a 
mean  very  close  to  the  point  X  with  only  the  appropriate 
bearings  from  the  other  stations  contributing  to  it.   Now 
if  the  computed  standard  deviation  is  small  enough  to  be 
less  than  a  location  confidence  limit,  X  is  called  a  pro- 
bable target  location.   The  observations  x,  ,  x2 ,  x^,  x^  are 
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flagged  as  correlated.   These  observations  are  deleted 
from  further  computations.   Then  the  next  observation  from 
station  I  is  selected  and  crossings  with  the  remaining 
bearing  lines  from  all  other  stations  are  calculated 
(Fig.  4c) .   Here  again,  by  computing  a  mean  location  and 
associated  standard  deviation,  and  rejecting  distant  cross- 
ings, the  computed  mean  should  converge  to  a  location  close 
to  Y  with  a  standard  deviation  smaller  than  the  preassigned 
limit.   Then  for  the  last  case,  ideally  only  bearings  lead- 
ing to  point  Z  are  left. 
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•Fig.  4.  Steps  for  target  location  by  "cluster"  points 
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III.   PRESENTATION  OF  RESULTS 

A.   TARGET  LOCATION 

For  the  computer  simulation,  two  target  configurations 
were  taken.   In  the  first,  three  targets  were  placed 
approximately  18  nmi  apart  on  a  line  parallel  to  the 
heading  of  the  observer  aircraft.   In  the  second,  six  tar- 
gets were  placed  on  a  circle  of  approximately  50  nmi 
diameter,  with  the  center  of  the  circle  approximately 
125  nmi  from  the  path  of  the  observer.   As  a  final 
aggravation,  the  probability  of  detection  of  any  target  was 
reduced  to  1/2;  i.e.,  only  three  targets  were  emitting  at 
any  instant.   A  random  sequence  determined  which  three  out 
of  the  six  targets  were  observed  at  each  station. 

In  the  simulation  the  "cluster"  method  was  first  tested, 
using  target  configuration  one;  Fig.  5  shows  the  results 
obtained.   Only  two  out  of  the  three  targets  were  found. 
The  position  errors  were  less  than  ten  miles.   The  third 
target  could  be  recovered  by  taking  a  second  look  at  the 
still  uncorrelated  observations.   To  get  an  estimate  on 
the  spread  of  the  computed  target  locations  that  could 
be  expected,  a  Monte  Carlo  simulation  was  done.   The  result 
is  shown  in  Fig.  6;  the  greatest  error  was  approximately 
50  nmi,  and  the  center  target  was  missed  twice.   The  com- 
puted mean  location  from  ten  runs  was  on  the  far  side  for 
all  three  targets  with  an  error  of  about  7  nmi. 
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Fig.  5.  Target  location  by  "cluster"  method,  single  run. 
Three  targets  are  approximately  18  NM  apart;  4  observation 
points  are  also  approximately  18nm  apart.  the  distance 
between  observations  and  targets  is  approximately  150  nh. 
Also  shown  are  the  noisy  observations  made  at  each  station 

+  true  target  location 

a  calculated  "cluster"  target  location 
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Fig,  6,  Ten  runs  of  target  location  by  "cluster"  method 

+  true  target  location 

a  computed  target  location/  individual  run 

b  mean  target  location  from  all  runs 

x  observation  points 
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Fig.  7.  Target  locations  obtained  from  "cluster"  method, 
Six  targets  are  placed  on  a  circle  of  25  NM  radius.  The 
center  of  the  circle  is  125  NH  from  the  base  line  of  obser- 
vation STATIONS. 

+  TRUE  TARGET  LOCATION 

A  COMPUTED  TARGET  LOCATION 

X   LOCATION  OF  OBSERVATION  POINTS 
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The  results  using  the  second  target  configuration  show 
only  one  of  the  computed  locations  is  reasonably  close  to 
the  corresponding  target,  and  that  looking  for  cluster 
points  is  not  very  effective  in  a  multi-emitter  environment, 
The  number  of  false  locations  becomes  so  high  that  the  pos- 
sibility of  random  clusters  formed  is  too  high  to  make  the 
occurrence  of  clusters  a  good  criterion  for  target  location. 
Therefore,  a  new  method  was  developed  using  a  Kalman  filter, 
and  a  correlation  procedure  to  obtain  refined  estimates  of 
the  observation  angles  from  the  measurements  taken. 
The  result  is  shown  in  Fig.  8;  all  positions  except  one  are 
closer  than  10  nmi  to  the  true  target  position.   To  allow 
a  comparison  between  the  two  methods ,  the  same  distance 
from  the  first  to  the  last  observation  point  was  chosen. 
In  Fig.  9  the  tracking  distance  is  doubled;  this  gave  a 
considerable  improvement.   The  largest  error  now  was  4.3  nmi 
with  3  positions  closer  than  2  nmi  to  the  true  target 
positions . 

As  a  final  test,  the  probability  of  detection  was 
reduced  to  1/2.   Fig.  10  shows  the  results  obtained;  given 
on  the  graph  are  the  initial  bearing  line  (smoothed)  where 
a  new  track  was  initialized;  the  estimates  of  bearing  lines 
at  the  "last  station;  the  target  location  based  on  only 
these  two  bearings  (0);  the  location  based  on  all  estimates 
available  (□);  the  true  target  locations  (+)  ;  and  the 
observation  points  (X) .   The  result  reflect  the  missing 
observations.   However,  the  computed  target  positions  are 
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Fig,  8,  Target  locations  obtained  when  using  Kalman  filter 
on  observation  angles.  the  bearing  lines  shown  are  the 
smoothed  estimates  obtained  at  the  first  and  last  observation 
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Fig,  9.  Target  location  using  Kalman  filter  with  20  obser- 
vation stations,  All  targets  are  observed  at  all  stations, 

+  true  target  location 

□  mean  target  location  using  all  available  bearing 
estimates 

0  target  location  using  first  and  last  bearing  estimate 
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Fig,  10.  .Target  location  using  Kalman  filter  with  20  obser- 
vation stations.  Only  3  targets  are  observed  at  any  one 
station.  The  graph  shows  the  estimate  of  the  initial  obser- 
vation AND  THE  LAST  OBSERVATION  OF  EACH  TRACK. 
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D  MEAN  TARGET  LOCATION  BASED  ON  ALL  AVAILABLE  ESTIMATES 

0   TARGET  LOCATION  BASED  ON  FIRST  AND  LAST  BEARING  ESTIMATE 
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still  close  enough  to  the  true  positions  to  be  useful  as 
apriori  target  locations. 

B.   TARGET  CORRELATION 

A  source  of  error  arises  in  the  target  correlation. 
Two  conflicting  requirements  make  a  compromise  necessary. 
First,  to  reduce  the  possibility  of  correlating  with  wrong 
observations,  it  is  desirable  to  select  a  small  gate  size, 
This,  however,  increases  the  possibility  of  missing  observa- 
tions which  should  be  correlated  with  a  track  but  happen 
to  be  noisier  than  expected.   This  is  undesirable,  because 
an  uncorrelated  observation  will  cause  initialization  of  a 
new  track  very  close  to  the  already  existing  track.   Now 
these  two  tracks  will  both  be  correlated  to  observations 
made  from  the  same  target,  and  these  observations  will  be 
split  by  the  two  tracks  —  degrading  the  quality  of  both. 
So,  false  initialization  of  new  tracks  is  even  less  desira- 
ble than  allowing  the  possibility  of  a  false  correlation. 
The  perfoirmance  of  the  correlation  method  was  evaluated 
for  several  test  runs.   Tables  I  and  II  show  two  typical 
cases.   The  numbering  of  tracks  and  observations  corresponds 
to  the  numbering  of  the  targets  shown  in  Fig.  7. 

The  percentage  of  wrong  correlations  of  20%  and  33.3% 
is  quite  high  and  degrades  the  filter  performance  seriously. 
Further  studies  should  be  done  to  find  a  more  effective 
way  of  correlation.   Even  when  noisefree  observations 
were  used,  the  noise  introduced  by  the  filter  due  to  the 
approximation  of  the  initial  angular  rate  and  the  closeness 
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TABLE  I 

Correlation  result  from  20  observation  points  with  all 
6  targets  observed  at  each  station.  all  tracks  are 
initialized  at  the  first  station. 
120  observations 
24  wrong  correlations  —  ►  20% 
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TABLE  II 

Correlation  result  from  20  observation  points  with  3 
targets  observed  at  each  station,  the  number  in  paren- 
thesis gives  the  station  at  which  the  track  was  initial- 
IZED, Track  numbers  are  assigned  equal  to  the  target 

NUMBER  THAT  HAS  CONTRIBUTED  THE  MOST  OBSERVATIONS  TO 
THIS  TRACK. 

60  OBSERVATIONS 

20  WRONG  CORRELATIONS *>       33.3% 
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of  the  targets  considered  was  sufficient  to  cause  correlation 

errors.   This  was  particularly  so  in  the  low  probability 

of  detection  case,  when  targets  were  emitting  close  together, 

and  only  one  track  was  initialized.   Then  this  track  tends 

to  pull  in  the  observations  from  both  targets,  unless 

they  are  observed  simultaneously  at  the  same  station. 

C.   STATISTICAL  EVALUATION 

Finally,  a  Monte  Carlo  simulation  was  performed  using  a 
single  target  (andideal  correlation)  to  verify  that  bias-free 
estimates  are  generated,  and  that  the  quality  of  the  esti- 
mates are  properly  represented  by  the  analytical  covariances 
obtained.   The  evaluation  is  done  for  the  Kalman  filter, 
with  and  without  the  correction  in  the  $(2,2)  term.   In 
both  cases,  the  covariance  sequences  of  the  estimates  of 
bearing  angle  and  angular  rate  are  reasonably  close  to 
the  analytical  values  (Fig.  11,  12);  but  without  the  correc- 
tion term,  a  bias  is  built  up  in  the  estimate  of  bearing 
angle  and  the  estimate  of  the  angular  rate.   With  the 
correction  term,  no  bias  is  present,  and  the  mean  error 
of  the  estimate  approaches  zero  (Fig.  2,3). 
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Fig.  11,  Variance  of  bearing  angle  estimate  obtained  from 
a  Monte  Carlo  simulation  on  the  Kalman  filter  with  a  single 
target  and  perfect  correlation.  First  observation  is  at 
time  zero. 
d  analytical  variance  of  the  angle  estimate 

+  experimental  variance  of  the  angle  estimate  when  using 
nonlinear  1  matrix  in  kalman  filter  prediction  equation 

x  experimental  variance  of  the  angle  estimate  when  using 
linear  1  matrix  in  kalman  filter  prediction  equation. 
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Fig.  12,  Variance  of  estimate  of  angular  rate  obtained  from 
a  Monte  Carlo  simulation  on  the  Kalman  filter  with  a  single 
target  and  perfect  correlation. 

□  analytical  variance  of  the  rate  estimate 

+  experimental  variance  of  the  rate  estimate  when  using 
nonlinear  *_  matrix 

x  experimental  variance  of  the  rate  estimate  when  using 

LINEAR  i  MATRIX 
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IV.  CONCLUSIONS 

The  results  from  the  computer  simulation  showed  that 
a  straight  forward  search  for  cluster  points  is  not  effec- 
tive in  a  high-density,  multiple-emitter  environment.   This 
is  so  becuase  if   m  targets  are  emitting,  and  there  are 
n  observations,  only   1/m  of  all  possible  locations  are 
true  locations—all  others  are  wrong!   This  ratio  of  true 
to  faulty  locations  can  be  reduced,  and  this  was  done  in  the 
computer  simulation  by  first  eliminating  locations  on  the 
wrong  side  of  the  flight  path,  eliminating  locations  too 
distant  to   be  considered  meaningful,  and  by  considering 
locations  along  selected  bearings  only.   The  high  number  of 
faulty  locations  still  remaining  is  sufficient  to  generate 
erroneous  clusters  to  which  the  process  may  converge  and 
lead  to  false  locations . 

VThen  a  Kalman  filter  was  used  to  process  the  measure- 
ments, the  critical  point  was  the  correlation  procedure 
which  selects  the  measurements  and  assigns  them  to  a  track. 
If  this  selection  was  done  perfectly,  the  filter  would  have 
only  the  proper  data  as  its  input  and  would  produce  an  opti- 
mal estimate  of  the  systems  state.   As  correlation  errors 
occur,  wrong  data  is  fed  into  the  filter  and  performance 
deteriorates.   This  indicates  that  the  overall  system  is 
only  as  good  as  the  correlation.   For  this  reason,  the 
correlation  procedure  is  quite  elaborate;  adaptive  gating 
was  introduced  and  five  rules  are  used  to  handle  the 
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multiple  correlation  problem  as  efficiently  as  possible. 
However,  it  must  be  pointed  out  that  the  two  requirements 
for  a  good  correlation  method — correlation  and  discrimination- 
in  a  sense  contradict  each  other  and  cannot  both  be  maximized 
at  the  same  time.   In  the  initial  correlation,  a  greater 
weight  is  put  on  the  correlation;  and  the  4  rejection  rules 
provide  discrimination,  but  errors  do  occur  during  this 
rejection  process,  and  the  high  error  percentage  obtained 
suggests  that  a  better  correlation  method  could  be  found 
which  then,  in  turn,  would  improve  the  filter  performance 
and  the  quality  of  locations  obtained. 


38 


APPENDIX  A 
KALMAN  FILTER  BLOCK  DIAGRAM  AND  INITIALIZATION 

Recursive  gain  and  covariance  equations  for  filtering: 


2K   =    PK/K-l^K[^-KPK/K-l^K+    V 


PT,  ,„    =    P 


hT,P, 


K/K   :        K/K-l         2.KiiK   K/K-l 


P  =    0  P  $  +0 

*K+1/K  K+1,K    K/K^K+1,K         UK 


Initial   values : 


1/0 


Hk 


"  4 
10    0 


0    1 


[1    0] 


K  ,K-1 


1    T 


0    1 


RK=  1 


QK  = 


0    0 


0    0 


for  stationary  targets. 
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Recursive  gain  and  covariance  equations  for  smoothing: 


T     -1 
n    =n     p       $     p 

Ul/K    Ul/K-1  K-1/K-1*K,K-1   K/K-l 


P,  ,„   =  P. 


h„P, 


1/K  "  Fl/K-1  "  U1/K?K-K  K/K-1U1/K 


Initial  values: 

1    0 


D 


1/1 


1/0 


0  1 

4 

10*  0 

0  1 


Initial  state,  values: 


1/1 


v 

—  sin  z, 

R      1 


v  =  velocity  of  aircraft 

R  =  150  NM  assumed  initial  distance  to  target 


Note:  In  the  Kalman  filter  equation  a  non-linear  $  matrix 
is  used ;when  calculating  gains  and  covariance s,  the  linear 
approximation  of  the  $  matrix  is  used. 
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Kalman  filter  gain  and  covariance  equations  in  scalar  form. 


*k(1)  -  Vi(1'1,/IVifl'1)  +  Rl 


V2>  ■  Vi(1'2)/IVill,1)  +  R1 


PK+1/K(1,D    =   gKd)R  +  gK(2)  [2R-PK/K_1(1,2)]    +  pk/r_1(2,2) 


PK+1/K(1'2)    "    gK(2,IR-PK/K-l(1'2)1    +    PK/K-1(2'2) 


PK+1/K(2'2)    =    PK/K-1(2'2)     "    PK/K-1(1'2)9k(2> 
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Fig,  13.  Block  diagram  of  Kalman  filter 
and  smoothing  solution 
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COMPUTER  RESULTS  FOR  SIX  TARGETS  AND  50%  DETECTION 
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BEARING    AMD    DISTANCE    FROM    STATION       1    TO    TARGET    INDICATED 


1 
2 
3 

% 

6 


305.23 
309.14 
303.82 
295.56 
290.31 
294.11 


126.36 
149.15 
170.16 
171.57 
151.83 
127.65 


BEARING  AND  DISTANCE  FROM  STATION   2  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


302.95 
307.31 
302.11 
293.73 
288.16 
291.60 


122.99 
145.43 
166.89 
169.07 
149.85 
125.31 


BEARING  AND  DISTANCE  FROM  STATION   3  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


300.55 
305.38 
300.33 
291.84 
285.95 
289.01 


119.83 
141.87 
163.78 
166.74 
148.09 
123.23 


BEARING  AND  DISTANCE  FROM  STATION   4  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


298.02 
303.36 
2  98.49 
289.90 
283.70 
286.33 


116.89 
138.49 
160.84 
164.60 
146.55 
121.41 


BEARING  AND  DISTANCE  FROM  STATION   5  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


295.36 
301.23 
296.57 
287.92 
281.40 
283.58 


114.20 
135.27 
158.06 
162.66 
145.25 
119.85 


BEARING  AND  DISTANCE  FROM  STATION   6  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


292.58 
299.01 
294.60 
285.88 
2  79.06 
280.76 


111.76 
132.27 
155.48 
160.91 
144.18 
118.59 
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BEARING  AMD  DISTANCE  FROM  STATION   7  TO  TARGET  INDICATED 

1  289,68  109.60 

2  296.69  129.47 

3  292.56  153.07 

4  283.81  159.38 

5  276.70  143.37 

6  277.89  117.63 


BEARING  AND  DISTANCE  FROM  STATION   8  TO  TARGET  INDICATED 

1  286.68  107.73 

2  294.27  126.88 

3  290.45  150.87 

4  281.70  158.06 

5  274.31  142.78 

6  274.98  116.94 


BEARING  AND  DISTANCE  FROM  STATION   9  TO  TARGET  INDICATED 

1  283.57  106.16 

2  291.75  124.54 

3  288.29  148.88 

4  279.56  156.95 

5  271.90  142.46 

6  272.04  116.58 


BEARING    AND    DISTANCE    FROM    STATION    10    TO    TARGET    INDICATED 

1  280.39  104.92 

2  289.15  122.44 

3  286.07  147.11 

4  277.38  156.06 

5  269.49  142.39 

6  269. C9  116.52 


BEARING  AND  DISTANCE  FROM  STATION  11  TO  TARGET  INDICATED 

1  277.14  104.00 

2  286.45  120.60 

3  283.80  145.56 

4  275.19  155.40 

5  267.08  142.57 

6  266.15  116.77 


BEARING    AND    DISTANCE    FROM    STATION    12    TO    TARGET    INDICATED 

1  273.84  103.42 

2  283.68  119.04 

3  281.49  144.25 

4  272.98  154.97 

5  264.68  143.00 

6  263.23  117.33 
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BEARING  AND  DISTANCE  FROM  STATION  13  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


270.51 
280.85 
279.14 
270.77 
262.29 
260.34 


103.20 
117.76 
143.17 
154.78 
143.68 
118.18 


BEARING  AND  DISTANCE  FROM  STATION  14  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


267.18 
277.96 
276.75 
268.55 
2  59.94 
257.50 


103.32 
116.79 
142.35 
154.82 
144.61 
119.35 


BEARING  ANO  DISTANCE  FROM  STATION  15  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


263.87 
275.02 
274.35 
266.33 
257.62 
254.72 


103.79 
116.11 
141.76 
155.09 
145.78 
120.78 


BEARING  AND  DISTANCE  FROM  STATION  1.6  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


260.60 
2  72.06 
271.92 
264.13 
255.34 
252.01 


104.60 
115.74 
141.43 
155.58 
147.17 
122.49 


BEARING  AND  DISTANCE  FROM  STATION  17  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


257.40 
269.09 
269.49 
261.94 
253.10 
249.39 


105.75 
115.68 
141.36 
156.32 
148.81 
124.48 


BEARING  AND  DISTANCE  FROM  STATION  18  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


254.27 
266.13 
267.07 
259.78 
250.92 
246.85 


107.21 
115.93 
141.54 
157.26 
150.66 
126.72 
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BEARING  AND  DISTANCE  FROM  STATION  19  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


251.23 
263.19 
264.65 
257.64 
248.80 
244.40 


108.99 
116.49 
141.97 
158.44 
152.73 
129.19 


BEARING  AND  DISTANCE  FROM  STATION  20  TO  TARGET  INDICATED 


1 
2 
3 
4 
5 
6 


248.30 
260.28 
262.25 
2  5  5.54 
246.73 
242.05 


111.06 
117.35 
142.66 
159.83 
155.00 
131.90 
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ESTIMATES  AT  STATION  1 

TRACK       THT  TH  K/K  TD  K/K      TH  1/K 

1        305.2349  306.4121  -1.8445  306.4121 

4  295.5591  2S5.1960  -2.0739  295.1960 

5  290.3093  290.5386  -2.1463  290.5386 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

304.567 

9999.500 

199.200 

3.500 

10.000 

293.122 

9999.500 

199.200 

3.500 

10.000 

288.392 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

293.828 

10000.000 

200.000 

3.200 

10.000 

290.897 

10000.000 

20C.000 

3.200 

10.000 

303.179 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
145000000 


4  0   10 

5  Oil 
3        10   0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1450000000 


4  0       10 

5  0      0       1 
3  10       0 


ESTIMATES  AT  STATION   2 

TRACK       THT         TH  K/K  TD  K/K      TH  1/K 

1        302.9536  303.6416  -2.3862  305.9487 

4  293.7261  293.5923  -1.9044  295.4312 

5  288.1587  290.0620  -1.3704  291.3738 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRE  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10. 000 

301.255 

9999.500 

1^9.200 

3.500 

10.000 

291.688 

9999.500 

199.200 

3.500 

10.000 

288.691 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

307.480 

10000.000 

20C.000 

3.200 

10.000 

285.740 

10000.000 

20C.000 

3.200 

10.000 

290.489 

INITIAL  CORRELATION  MATRIX 

TRACKS 


OBSERVATIONS 


0   0   0   0   0   0 


2 

0 

0 

0 

5 

0 

0 

1 

6 

0 

1 

1 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1450000000 


2  0       0       0 

5  0      0       1 

6  0      10 


ESTIMATES  AT  STATION   3 

TRACK       THT         TH  K/K  TD  K/K      TH  1/K 

1  300.5466    301.2554  -2.5068  305.9487 

4  291.8401    290.8882  -2.3523  295.4307 

5  285.9539    286.7239  -2.3736  291.3733 

2  305.3792    307.4797  -1.9072  307.4797 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 

9999.500  199.200 

9999.500  199.200 

9999.500  199.200 

10000.000  20C.000 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


3.500 

9.000 

298.749 

3.500 

10.000 

288.536 

3.500 

10.000 

284.350 

3.200 

10.000 

305.573 

10000.000 

200.000 

3.200 

10.000 

284.413 

10000.000 

200.000 

3.200 

10.000 

286.115 

10000.000 

200.000 

3.200 

10.000 

298.096 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


6  0       110 

5  0      110 

1  10      0      0 

FINAL  CORRELATION  MATRIX 


OBSERVATICNS 


TRACKS 
1452000000 


6 

0 

0 

1 

0 

5 

0 

1 

0 

0 

1 

1 

0 

0 

0 

ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 

298.0159 
289.9031 
283.6997 
303.3562 

298.3137 
287.0227 
264.3892 
305.5725 

-2.8423 
-3.0166 
-2.4085 
-1.9980 

305.9485 
295. 7332 
291.3655 
307.4797 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 

9999.500  199.200 

9999.500  199.200 

9999.500  199.200 

10000.000  20C.000 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


3.500 

10.000 

295.471 

3.500 

10.000 

284.006 

3.500 

10.000 

281.980 

3.200 

9.000 

303.574 

10000.000 

200.000 

3.200 

10.000 

292.960 

10000.000 

200.000 

3.200 

10.000 

299.960 

10000.000 

20C.000 

3.200 

10.000 

282.852 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


1 

1 

0 

0 

0 

2 

0 

0 

0 

1 

6 

0 

1 

1 

0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


1  10   0   0 

2  0   0   0   1 
6        0   0   10 


ESTIMATES  AT  STATION   5 

TRACK       THT         TH  K/K  TD  K/K  TH  1/K 

1  295.3552  293.9016  -3.5951  306.2622 

4  287.9155  284.0059  -3.0959  295.7332 

5  281.4001  282.4714  -2.2949  291.2229 

2  301.2344  301.1648  -3.2866  306.2742 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRE  PW        COUNT     DOA 


9999.500 

1^9.200 

3.500 

10.000 

290.306 

9999.500 

199.200 

3.500 

9.000 

280.910 

9999.500 

199.200 

3.500 

10.000 

280.177 

10000.000 

200.000 

3.200 

10.000 

297.878 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

200.000 

3.200 

10.000 

277.730 

10000.000 

200.000 

3.200 

10.000 

283.602 

10000.000 

200.000 

3.200 

10.000 

292.841 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


5  0      110 

4  0      110 

1  10       0      0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


5 

0 

0 

1 

0 

4 

0 

1 

0 

0 

1 

1 

0 

0 

0 

ESTIMATES  AT  STATION   6 

TRACK  THT         TH  K/K  TD  K/K      TH  1/K 

1  292.5771  291.7351  -3.3142  305.8477 

4  285.8843  282.4272  -2.6802  295.2927 

5  279.0637  278.9421  -2.6546  291.6421 

2  299.0105  297.8782  -3.4860  306.2742 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

288.421 

9999.500 

199.200 

3.500 

10.000 

279.747 

9999.500 

199.200 

3.500 

10.000 

276.287 

0000.000 

20C.0C0 

3.200 

9.000 

294.392 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

290.938 

10000.000 

200.000 

3.200 

10.000 

283.122 

10000.000 

20C.000 

3.200 

10.000 

296.286 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


1  10   0   1 
4        0   10   0 

2  0   0   0   1 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


1  10       0      0 
4  0       10      0 

2  0      0      0       1 


ESTIMATES  AT  STATION   7 

TRACK       THT         TH  K/K  TD  K/K  TH  1/K 

1  289.6799    289.6907  -3.1110  305.4155 

4  283.8108    281.4497  -2.2746  294.7136 

5  276.6963    276.2874  -2.6817  291.6421 

2  296.6895    295.6548  -3.0584  306.2739 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PVi  COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

286.580 

9999.500 

199.200 

3.500 

10.000 

279.175 

9999.500 

199.200 

3.500 

9.000 

273.605 

10000.000 

20C.000 

3.200 

10.000 

292.596 

INPUT  OBSERVATIONS 

FREQU       PRF  PW  COUNT  DOA 

10000.000    200.000      3.200  10.000  274.281 

10000.000    20C.000      3.200  10.000  275.103 

10000.000    20C.000      3.200  10.000  288.125 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


5  0   0   10 

6  0   0   10 
1        10   0   0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452000000 


5  0   0   10 

6  0   0   0   0 
1        10   0   0 


ESTIMATES  AT  STATION   8 

TRACK       THT         TH  K/K  TD  K/K  TH  1/K 

1  286.6750  287.2795  -3.0562  305.1560 

4  281.6995  279.1750  -2.3038  294.7136 

5  274.3062  273.9114  -2.6289  291.5286 

2  294.2708  292.5962  -3.1943  306.2739 

6  274.9802  275.1030  -2.2990  275.1030 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

284.223 

9999.500 

199.200 

3.500 

9.000 

276.871 

9999.500 

199.200 

3.500 

10.000 

271.282 

10000.000 

20C.000 

3.200 

9.000 

289.402 

10000.000 

20C.000 

3.200 

10.000 

272.804 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 

10000.000    20C.000 
10000.000    20C.000 


10000.000    20C.000 
INITIAL  CORRELATION  MATRIX 


3.200 

10.000 

285.256 

3.200 

10.000 

290.960 

3.200 

10.000 

272.636 

OBSERVATIONS 

1 
2 

6        0   0   10   1 

FINAL  CORRELATION  MATRIX 


TRACKS 

1 

4   5   2 

6 

1 

0 

0   0   1 
0   0   1 

0 
0 

0   0   0   0   0 


OBSERVATIONS 


TRACKS 
1452600000 


1  10   0   0   0 

2  0   0   0   10 
6        0   0   0   0   1 


ESTIMATES  AT  STATION  9 

TRACK  THT  TH  K/K  TD  K/K  TH  1/K 

1  283.5730  284.6465  -3.0573  304.9905 

4  279.5559  276.8711  -2.3261  294.7136 

5  271.9006  271.2822  -2.6343  291.5286 

2  291.7529  290.3755  -2.9371  306.0789 

6  272.0425  272.6921  -2.3635  275.0469 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10,000 

281.589 

9999.500 

199.200 

3.500 

8.000 

274.545 

9999.500 

1°9.200 

3.500 

9.000 

268.648 

10000.000 

20C.000 

3.200 

10.000 

287.438 

10000.000 

2CC.000 

3.200 

10.000 

270.329 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

280.466 

10000.000 

20C.000 

3.200 

10.000 

289.489 

10000.000 

20C.000 

3.200 

10.000 

287.626 

INITIAL  CCRRELATION  MATRIX 


TRACKS 

ATICNS 

1 

4 

5 

2 

6 

1 
2 
3 

1 
0 
0 

0 
0 

0 

0 
0 
0 

0 

1 
1 

0 
0 
0 

0   0   0   0   0 


FINAL  CORRELATION  MATRIX 


OBSERVATIONS 

1 
2 
3  0      0       0      10 


TRACKS 
1452600000 


1  10   0   0   0 

2  0   0   0   0   0 


ESTIMATES  AT  STATION  10 

TRACK       THT         TH  K/K  TD  K/K  TH  1/K 

1  280.3870  281.1697  -3.1954  305.1599 

4  277.3848  274.5449  -2.3411  294.7136 

5  269.4885  268.6477  -2.6286  291.5286 

2  289.1460  287.5442  -2.9981  306.0481 

6  269.0938  270.3286  -2.3646  275.0469 
2  289.1460  289.4893  -2.2181  289.4893 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRE  PW        COUNT     DOA 

9999.500  199.200 

9999.500  199.200 

9999.500  199.200 

10000.000  20C.000 

10000.000  2C0.O0O 

10000.000  20C.000 

INPUT  OBSERVATIONS 

FREOU       PRF  PW        COUNT     DOA 


3.500 

10.000 

277.974 

3.500 

7.000 

272.204 

3.500 

8.000 

266.019 

3.200 

10.000 

284.546 

3.200 

9.000 

267.964 

3.200 

10.000 

287.271 

10000.000 

20C.000 

3.200 

10.000 

265.908 

10000.000 

2C0.000 

3.200 

10.000 

277.801 

10000.000 

20C.000 

3.200 

10.000 

267.521 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


5  0      0       10  10 
1                    10       0      0  0      0 

6  0      0       10  10 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


5  0   0   10   0   0 
1        10   0   0   0   0 

6  0   0   0   0   10 


ESTIMATES  AT  STATION  11 


ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 
6 
2 

277.1357 
275.1921 
267.0781 
286.4541 
266.1504 
286.4541 

277.9148 
272.2036 
265.9734 
284.5459 
267.6685 
287.2710 

-3.2543 
-2.3485 
-2.6207 
-3.0795 
-2.5043 
-2.2715 

305.1843 
294.7136 
291.5461 
306.0481 
275.0466 
289.4893 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

274.660 

9999.500 

199.200 

3.500 

6.000 

269.855 

9999.500 

199.200 

3.500 

10.000 

263.353 

10000.000 

200.000 

3.200 

9.000 

281.466 

10000.000 

200.000 

3.200 

10.000 

265. 164 

10000.000 

20C.000 

3.200 

9.000 

284.999 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

200.000 

3.200 

10.000 

265.505 

10000.000 

200.000 

3.200 

10.000 

274.497 

10000.000 

20C.000 

3.200 

10.000 

280.377 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


5  0  0  10  10 
1  10  0  0  0  0 
3        0   0   0   10   1 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
145262000 


5  0  0  0  0  10 
1  10  0  0  0  0 
3        0   0   0   10   0 


ESTIMATES    AT    STATION    12 
TRACK 


1 
4 
5 
2 
6 
2 


THT 

TH 

K/K 

TD  K/K 

TH 

1/K 

273.8369 

274, 

,6086 

-3.2929 

305, 

,2207 

272.9841 

269. 

.8550 

-2.3480 

294, 

.7136 

264.6780 

263. 

,3525 

-2.5927 

291. 

,5461 

283.6846 

280, 

.9167 

-3.2905 

306. 

,2344 

263.2268 

265. 

.3770 

-2.4014 

275, 

,0039 

283.6846 

284, 

,9993 

-2.3198 

289, 

,4893 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 

9999.500  199.200 

9999.500  199.200 

9999.500  199.200 

10000.000  20C.000 

10000.000  200.000 

10000.000  200.000 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


3.500 

10.000 

271.316 

3.500 

5.000 

267.507 

3.500 

9.000 

260.760 

3.200 

10.000 

277.626 

3.200 

10.000 

262.975 

3.200 

8.000 

282.679 

10000.000 

20C.000 

3.200 

10.000 

278.386 

10000.000 

200.000 

3.200 

10.000 

259.362 

10000.000 

20C.000 

3.200 

10.000 

270.416 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


3  0      0       0       10       1 
6                    0      0       10       10 

4  110      0       0      0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


3  0   0   0   10   0 
6        0   0   10   0   0 

4  10   0   0   0   0 


ESTIMATES  AT  STATION  13 

TRACK  THT         TH  K/K  TD  K/K  TH  1/K 

1  270.5105  271.0515  -3.3342  305.3345 

4  270.7659  267.5068  -2.3396  294.7136 

5  262.2949  260.2375  -2.6440  291.7573 

2  280.8472  277.9705  -3.2648  306.1067 

6  260.3367  262.9753  -2.3766  275.0039 
2  280.8472  282.6794  -2.3621  289.4893 
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ORIGINAL  MATRIX  X  OF  PREDICTEO  TRACKS 

FREQU       PRF  PW        COUNT     DOA 

9999.500  199.200 

9999.500  199.200 

9999.500  199.200 

10000.000  2CC.000 

10000.000  20C.000 

10000.000  20C.000 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


3.500 

10.000 

267.717 

3.500 

4.000 

2  6  5.167 

3.500 

10.000 

257.594 

3.200 

10.000 

274.706 

3.200 

9.000 

260.599 

3.200 

7.000 

280.317 

10000.000 

20C.000 

3.200 

10.000 

266.405 

10000.000 

2C0.OOO 

3.200 

10.000 

276.800 

10000.000 

20C.000 

3.200 

10.000 

257.448 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1  110      0       0      0 

3  0      0       0       10       1 

6  0      0       10       10 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
145262000 


1  0      10      0      0      0 

3  0      0       0      10      0 

6  0      0      10      0      0 


ESTIMATES    AT    STATION    14 


ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 
6 
2 

267.1829 
268.5466 
259.9407 
277.9558 
257.4980 
277.9558 

267.7173 
265.7283 
257.5435 
275.5640 
260.5986 
280.3174 

-3.3187 
-2.1972 
-2.5980 
-3.1306 
-2.3439 
-2.3975 

305.3345 
294.5056 
291.7778 
305.7715 
275.0039 
289.4893 

61 


ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

9.000 

264.398 

9999.500 

199.200 

3.500 

10.000 

263.531 

9999.500 

199.200 

3.500 

10.000 

254.945 

10000.000 

20C.000 

3.200 

10.000 

272.433 

10000.000 

200.000 

3.200 

8.000 

258.255 

10000.000 

2CC.000 

3.200 

6.000 

277.920 

INPUT  OBSERVATIONS 

FREQU       PRF          PW  COUNT  DOA 

10000.000    20C.000      3.200  10.000  273.934 

10000.000    200.000      3.200  10.000  267.236 

10000.000    20C.000      3.200  10.000  257.372 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


2        0   0   0   10   1 

4  10   0   0   0   0 

5  0   0   10   10 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


2        0   0   0   10   0 

4  10   0   0   0   0 

5  0   0   0   0   10 


ESTIMATES  AT  STATION  15 

TRACK  THT         TH  K/K  TD  K/K  TH  1/K 

1  263.8740  265.1750  -3.1932  304.9941 

4  266.3318  263.5310  -2.1781  294.5056 

5  257.6194  254.9455  -2.5346  291.7778 

2  275.0237  272.9934  -3.0500  305.5442 

6  254.7192  257.7568  -2.4628  275.1479 
2  275.0237  277.9197  -2.4254  289.4893 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

261.982 

9999.500 

199.200 

3.500 

9.000 

261.353 

9999.500 

1^9.200 

3.500 

9.000 

252.411 

10000.000 

20C.000 

3.200 

10.000 

269.943 

10000.000 

200.000 

3.200 

10.000 

255. 2^4 

10000.000 

20C.CO0 

3.200 

5.000 

275.494 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

2  0C.000 

3.200 

10.000 

262.133 

10000.000 

20C.000 

3.200 

10.000 

256.917 

10000.000 

20C.000 

3.200 

10.000 

272.597 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1  110   0   0   0 
5        0   0   0   0   10 

2  0   0   0   10   1 

FINAL    CORRELATION    MATRIX 


OBSERVATIONS 


TRACKS 
145262000 


1  10   0   0   0   0 
5        0   0   0   0   10 

2  0   0   0   10   0 


ESTIMATES  AT  STATION  16 

TRACK  THT         TH  K/K  TD  K/K      TH  1/K 

1  260.60^0  262.0203  -3.1390  304.9768 

4  264.1274  261.3528  -2.1529  294.5056 

5  255.3379  252.4109  -2.4635  291.7778 

2  272.0647  270.8521  -2.9119  305.1655 

6  252.0108  256.1128  -2.1910  274.8694 
2  272.0647  275.4941  -2.4452  289.4893 
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ORIGINAL 

MATRIX  X  OF 

PREDICTED 

TRACKS 

FREQU 

PRF 

PW 

COUNT 

DOA 

9999.500 

1°9.200 

3.500 

10.000 

258. 

,881 

9999.500 

199.200 

3.500 

8.000 

259. 

,200 

9999.500 

199.200 

3.500 

8.000 

249, 

,947 

10000. 000 

200,000 

3.200 

10.000 

267. 

,940 

10000.000 

2CC.000 

3.200 

10.000 

253. 

,922 

10000.000 

2CC.000 

3.200 

4.000 

273, 

,049 

INPUT  OBSERVATIONS 
FREQU       PRF 


PW 


COUNT 


INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 

4 
5 
6 


TRACKS 
14       5       2       6       2 


110  0  0  0 
0  0  10  10 
0      0       10       10 


FINAL  CORRELATION  MATRIX 


DOA 


10000.000 

200.000 

3.200 

10.000 

261.553 

10000.000 

200.000 

3.200 

10.000 

251.953 

10000.000 

20C.000 

3. 200 

10.000 

250.358 

0   0   0   0 


OBSERVATIONS 

4 
5 
6 


TRACKS 
145262000 


0  10  0  0  0 
0  0  0  0  10 
0   0   10   0   0 


ESTIMATES  AT  STATION  17 
TRACK       THT         TH  K/K 


TD  K/K 


TH  1/K 


1 

257.3982 

258.8811 

-3.0714 

304.9768 

4 

261.9414 

260.1643 

-1.9332 

294.1287 

5 

253.1047 

250.0772 

-2.3646 

291.6860 

2 

269.0950 

267.9402 

-2.9012 

305.1655 

6 

249.3858 

253.0294 

-2.3437 

275.1997 

2 

269.0950 

273.0488 

-2.4563 

289.4893 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

9.000 

255. 810 

9999.500 

199.200 

3.500 

10.000 

258.231 

9999.500 

199.200 

3.500 

10.000 

247.713 

10000.000 

200.000 

3.200 

9.000 

265.039 

10000.000 

200.000 

3.200 

10.000 

250.686 

10000.000 

20C.000 

3.200 

3.000 

270.593 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

255.230 

10000.000 

200.000 

3.200 

10.000 

268. 186 

10000.000 

2CC.000 

3.200 

10.000 

260.174 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1        110   0   0  0 

3  0   0   0   10  1 

4  0   10   0   0  0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1        10   0   0   0   0 

3  0   0   0   0   0   1 

4  0   10   0   0   0 


ESTIMATES  AT  STATION  18 


ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 
6 
2 

254.2676 
259.7776 
250.9225 
266.1299 
246.8473 
266.1299 

255.6700 
258.9563 
247.7126 
265.0388 
250.6858 
268.9878 

-3.0017 
-1.7804 
-2.2846 
-2.8757 
-2.2765 
-3.2545 

305.0388 
293.8345 
291.6860 
305. 1655 
275.1997 
288.6863 
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ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

252.668 

9999.500 

199.200 

3.500 

1C.000 

257.176 

9999.500 

199.200 

3.500 

9.000 

245.428 

10000.000 

200.000 

3.200 

8.000 

262.163 

10000.000 

200.000 

3.200 

9.000 

248.409 

10000.000 

20C.000 

3.200 

10.000 

265.733 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

2CC.000 

3.200 

10.000 

251.941 

10000.000 

200.000 

3.200 

10.000 

263.135 

10000.000 

20C.000 

3.200 

10.000 

245.247 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 

1 
2 
6 

FINAL  CORRELATION  MATRIX 


TRA 

CKS 

1 

4 

5 

2 

6 

2 

1 

0 
0 

0 
0 
0 

0 
0 
1 

0 
1 
0 

1 

0 

1 

0 
1 
0 

0   0   0   0 


OBSERVATIONS 


TRACKS 
145262000 


1  10   0   0   0   0 

2  0   0   0   10   0 
6        0   0   10   0   0 


ESTIMATES  AT  STATION  19 


ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 
6 
2 

251.2313 
257.6436 
248.7965 
263.1853 
244.4004 
263.1853 

252.5027 
257.1758 
245.3749 
262.4705 
248.4093 
265.7332 

-2.9186 
-1.7552 
-2.2080 
-2.7851 
-2.204O 
-3.2269 

305.1133 
293.8345 
291.7087 
304.9487 
275.1997 
288.6863 

66 


ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS 

FREQU       PRF  PW        COUNT     DOA 


9999.500 

199.200 

3.500 

10.000 

249.584 

9999.500 

199.200 

3.500 

9.000 

255.421 

9999.500 

199.200 

3.500 

10.000 

243.167 

10000.000 

20C.000 

3.200 

10.000 

259.685 

10000.000 

2C0.OOO 

3.200 

8.000 

246.204 

10000.000 

200.000 

3.200 

9.000 

262.506 

INPUT  OBSERVATIONS 

FREQU       PRF  PW        COUNT     DOA 


10000.000 

20C.000 

3.200 

10.000 

247.938 

10000.000 

200.000 

3.200 

10.000 

263.600 

10000.000 

20C.000 

3.200 

10.000 

240.869 

INITIAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1        10   0   0  10 

3        0   0   0   0  0   1 

6        0   0   10  0   0 

FINAL  CORRELATION  MATRIX 


OBSERVATIONS 


TRACKS 
1452620000 


1  10  0  0  0  0 
3  0  0  0  0  0  1 
6        0   0   10   0   0 


ESTIMATES  AT  STATION  20 


ACK 

THT 

TH  K/K 

TD  K/K 

TH  1/K 

1 
4 
5 
2 
6 
2 

248.3002 
255.5435 
246.7299 
260.2766 
242.0506 
260.2766 

249.2296 
255.4206 
242.5381 
259.6853 
246.2044 
263.2356 

-2.8380 
-1.7272 
-2.1949 
-2.7359 
-2.1300 
-2.8191 

305.2737 

293.8345 
291.9839 
304. 9487 
275.1997 
288.6863 
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COMPUTER  PROGRAM  FOR  APRIORI  TARGET  LOCATION  AND  CORRELATION 


C   LIST  OF  VARIABLES 

C      THT(I,J)       TRUE  BEARING  ANGLE 

C      THETA(I,J)     MEASURED  BEARING  ANGLE 

C      TH1KK(I,J)     FILTERED  ESTIMATE  OF  BEARING  ANGLE 

C      TD1KK(I,J)     FILTERED  ESTIMATE  OF  ANGULAR  RATE 

C      THKJ)         SMOOTHED  ESTIMATE  OF  FIRST  BEARING  ANGLE 

C      E1(I,J)         ERROR  OF  TH1KK(I,J) 

C      NJ(I,J)        RANDOM  SEQUENCE  OF  TARGETS  DETECTED 

C 

C 

C   LIST  OF  SUBSCRIPTS 

C 

C      I  l.GE.I.LE.NSTA       STATION  INDEX 

C      K  1.GE.K.LE.20         K.F.  INDEX 

C   NOTE:   I  AND  K  WILL  IN  GENERAL  BE  THE  SAME 

C      J  l.GE. J.LE.NTAR       TARGET  INDEX 

C 

REAL*8  ITL( 12) ,LBL/'  '/ 

REAL  LABEL/1     •/ 

COMMON  /A/  NTR,NJ,TH1,TD1KK,NFRST  /B/TL AT , TLOT, I TL, X, Y 
DIMENSION  THT(  20,20) ,THET A( 20 , 20 ) , STEESTt 20,20) ,XX( 20) 
1TD1KM  20,20),GK(20)  ,GD(20) ,PKK( 20) ,NFRST(20) ,NJ(20,20) 
1E1(20,20),TH(20) ,TLA(20) ,TLO(20) , SLA(20) , SLC ( 20 ) , X ( 20) 
IT L AT (20) ,0BS<5,20) ,TRK(5,20) , GATE ( 5 ) , ADDGT ( 20 ) , YY ( 20) , 
1TLOT(20) ,NK1(20)  ,NTR(20) ,Y(20)  ,TH1KK(20,20)  ,TH1(20)  , 
1IC(20),SLA1(20),SL01(20) 
C 

DATA  GK  /l. , .6667, .6667, .625, .5637, .5045, .4532, .4099, 
X.3734, .342 5,. 3 161, .29  34, .2737,-2  563, . 24 1 1 , . 2275 , . 21 53 , 
X. 2044, .1945, .1855/ 

DATA  GD/O. ,. 33 32, .3 33 3,. 25,. 1818,. 13 52,. 103 5,. 08 141, 
X.C65  58, .0539,.  053  9, .03  32, .0328, .02  845 ,. 02  49 , .022, 
X. 01956,. 01 7 5,. 01 576,. 01 427/ 

DATA  PKK/1.414, 1.414, 1.29,1. 13, 1.007, .91,. 833,. 772, 
X. 724,. 6 8,. 64 5, . 61 5 , . 5 865 , . 564, . 542 , . 524 , . 507 , .497 , .47  7 
X.4645/ 
DATA  PI/3.141593  /, RD/. 0174533/ , NK1 /20*0/ , R/3959. / 
CATA  GATE/1. ,1 .,3.,  10.  ,2./ 
C 

IX=38159 

TDM=2.292 

READ  (5,102)  ( ITL( I) ,1=1,12) 

DO  1000  J=1T6 

1000  READ  (5,103)  ( OBS( I , J ) , TRK ( I , J ) , I =1 ,4 ) 
DO  J001  K=l,20 

1001  READ  (5,105)  ( NJ ( K , J  )  , J  =  l ,3 ) 


C 


READ  (5,105)  (NTR( J)  ,J=1,20) 

CALL  BEARG1  (THT, STEEST ,NSTA, NO , SLA, SLO ) 


C 

C    INITIALIZATION 


C 


1  =  1 

WRITE  (6,101) 
WRITE  (6,106)  I 
N0=3 

DO  10  J=1,NC 
J1=NJ( 1, J) 
NTR( J)=J1 
NFRST( J)=l 

CALL  GAUSS ( IX,STEEST(1,J),THT(1,J1) ,THETA( 1,J  )  ) 
TH1KK( 1 , J)=THETA( 1, J) 
TDIKMI,  J  J  =TDM*SIN(THETA(  1,  J)*RD) 
THK  J)=THETA(1  ,J) 
E1(1,J)=TH1KK( 1, J)-THT( 1, Jl) 

WRITE  (6,100)NTR( J ) ,THT(1,J1) ,TH1KK( 1 , J ) , TD IKK ( 1 , J ) 
STHK  J) 
10   CONTINUE 


68 


c 

C    LINEAR  KALMAN  FILTER  WITH  PREDET.  GAINSCHEDULE 


C 
C 

c 


NT1=N0 

DO  12  I=2,NSTA 

DO  11  J=1,NT1 


IF(J.GT.NO)  GOTO  1 

J1=NJ( I, J) 

CALL  GAUSS ( IX,STEEST(I,J) ,THT( I,J1) ,THETA( I, J) ) 

OBS(5, J)=THETA( I , J) 

I  K=I-NKMJ) 
C 

C    COMPUTE  PREDICTED  VALUE  FOR  CORRELATION 
C 

TRK(5, J)=TH1KK< I -1 , J ) +TD1KK ( I-1,J) 

ADDGT( J)=PKK(K)*2. 

II  CONTINUE 
WRITE  (6,101) 

CALL  CORREL  (I C,N0,NT1 ,OBS ,TRK,GATEt ADDGT ,NK1, I ) 

WRITE  (6t106)  I 

DO  12  J=1,NT1 
C 

C   COMPUTE  ESTIMATE 
C 

K=I-NK1 < J) 

E=0BS(5t J)-TRK(5,J) 

ET  =  0. 

ED  =  0. 

IF(  IC( J ) .EQ.O)  GOTO  110 

ET=GK(K)*E 

ED=GD(K )*E 
110  CONTINUE 

TH1KK( I , J)=TRK(5,J)+ET 

EPS=2.*TD1KK( I-1TJ)*RD/TAN(TH1KK( I, J)*RD) 

TDKKM=TD1KK(I-1,J) Ml.  +  EPS) 

TD1KK( I , J)=TDKKM+ED 

NTRJ=NTR(J) 

El ( I  ,J)=TH1KK(  I , J)-THT(I ,NTRJ) 

THK  J)=TH1(  J)+ET-(K-1)*ED 
C 

C         WRITE    RESULTS 
C 

WRITE(6t100)NTRJ,THT(I ,NTRJ)  ,TH1KK(  I , J)  , TD  IKK ( I , J ) 
GTHK  J) 
12       CONTINUE 
C 
C         COMPUTE    APRIORI    TARGET    LOCATION 


C 


WRITE(6,104) 

DO    121    1=1, NT1 

IF(TRK(4,I ) .LT.O )    GOTO    121 

TH(1)=TH1( I ) 

NFI=NFRST( I ) 

SLAK1  )  =  SLA(NFI  ) 

SLOK  1)=SLG(NF  I) 

NS=NSTA-NFI+1 

DO    120    J=2,NS 

K=NFI+J-1 

TH( J)=TH1KK(K,  I  ) 

SLAK J)=SLA(K) 

SLOK  J)=SLO(K) 

120  CONTINUE 

CALL  FIND  (TLA(  I),TLO( I),SIG,NS,TH,SLA1,SL01) 

NTRJ=NTR(I ) 

DELLA=TLA( I )-TLAT(NTRJ) 

DELLO=TLO( I )-TLOT(NTRJ) 

RNGERR=SQRT( (DELLA*60. ) **2+( DELL 0*50. )**2) 

WRITE  (6,107)  TLA( I),TLO( I) , DELL A T DELLO , RNGERR 

121  CONTINUE 


69 


C   PLOT  NOISY  BEARING  LINES 
C 

DO  1002  L=l,2 

DO  1002  1=1, NT1 

IF(TRK(4,I ) .LT.O)  GOTO  1002 

K=NFRST( I ) 

IFCL.EQ.2)  K=NSTA 

XX(1)=X(K) 

YY(1)=Y(K) 

THR=TH1 ( I) 

IFU.EQ.2)  THR  =  TH1KK(NSTA,  I  ) 

IF(THR.GT.180)  THR=360.-THR 

CA=ABS(TLO(  I)-SLO(K) )*RD 

SC=(90.-SLA(K) )*RD 

SB=(90.-TLA( I ) )*RD 

SB1  =  ARC0S(C0S( SB)*COS( SC )  +  SIN ( SB ) *SIN ( SC ) *COS  (  CA)  ) 

DO  1010  J=l,6 

SA=ARCOS(COS(SB)*COS( SC)+SIN(S3)*SIN< SC) *COS(THR*RD ) ) 
CB=ARCOS( (COS(SB)-COS(SA)*COS(SC) )/ (SIN( SA )*SIN( SC) ) ) 

IF(SLO(K).LT.TLO(I) )  GOTO  1008 

THA=CB-SLO(K)*RD 

GO  TO  1009 
10C8  THA=-CB-SLO(K)*RD 
10C9  CONTINUE 

THA=-THA/RD 

PH=90.-SA/RP 

CALL  SHIFT  (R , THA , PH, XX ( J+l )  , YY ( J+l ) ) 
1010  CONTINUE 

CALL  DRAW  ( 7  ,  XX , YY, 2, 0, LABE L , I TL , . 2 , . 2 ,5 , 5 , 2 , 2,9 , 10 ,0 , 
1002  CONTINUE 
C 
C    PLOT  COMPUTED  TARGET  POSITION 


C 


K=0 

DO  13  1=1, NT1 
IF(TRK(4,I ) .LT.O)  GOTO  13 
K  =  K+1 

CALL  SHIFT  ( 1 . , TLO(  I ) , TL AC  I ) , XX ( K ) , YY ( K ) ) 
13   CONTINUE 

CALL  DRAW  (  K ,XX , YY , 3,3 , LBL , ITL , . 1 , . 1 , 5, 5 , 2 , 2, 9, 10 ,0 , L ) 
STOP 

100  FORMAT  ( 117, 4X ,10F11.4) 

101  FORMAT  (1H1) 

102  FORMAT  (6A8) 

103  FORMAT  (8F10.5) 

104  FORMAT  (  1H1 , // / / 1 4X  ,  • TARLAT • , 4X , ' TARLON • , 4X , ■ DEL . LAT • , 
£3X,'DEL.L0N* ,3 X, • RNG. ERR ' , / ) 

105  FORMAT  (4012) 

106  FORMAT     ( /// 14X ,  « EST IMATE S    AT    STAT  I  ON «  , I  3 , //14X , • TRACK • 
6,6X,,THT,,8X,,TH    K/K,,5X,'TD    K/KS5X,«TH     1/K',/) 

1C7  FORMAT  (11X,8F10.5) 
END 

SUBROUTINE  SHI  FT ( R , TH ,PH, X, Y ) 

CATA  IFLAG  /l/ 

IF(  IFLAG. EQ.O)  GOTO  100 

IFLAG=0 

XO=TH 

YO  =  PH 

X=0. 

Y=.5 

RETURN 

100  CONTINUE 
X=TH-XO 
Y=PH-Y0+.5 
RETURN 
END 
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SUBROUTINE  BEARG1  ( THT , STEEST ,NST , NTAR , SL A, SLO ) 
C 

C   THIS  PROGRAM  GENERATES  BEARINGS  FROM  A  SET  OF  STATION 
C   AND  TARGET  POINTS.  THE  BEARINGS  ARE  GIVEN  AS  ANGLES  FROM 
C   0  TO  360  DEGREES,  WHERE  0  DEGREES  IS  TRUE  NORTH. 
C 

REAL*8  ITK12)  ,LBL/»         •/ 

REAL  LB  /•     «/ 

COMMON  /B/  TLA, TLO, ITL,XX,YY 

DIMENSION  TLA(2O),TLO(20),SLA(20),SLO(2O) , STDEV(20) , 
1 DIST (20, 20) T STEEST (20 ,20) , HD ( 20 ) , XX ( 2  0) , YY( 20 ) , X ( 20 ) , 
1THT(20,20) ,Y(20) 
C 

TW0PI=2.*3. 14159265 

PI=3. 141592 

RADIUS=10800./PI 

RADDEG=PI/180.0 
C 

C   READ  TARGET  POINTS 
C 

READ  (5,100)  NTAR 

DO  1  1  =  1, NTAR 

1  READ(5,102)  TLA(  I)  , TLO ( I ) 
C 

C   READ  STATION  POINTS 
C 

READ(5,100)  NST 

DO  2  1  =  1, NST 

2  READ(5,102)  SL A (  I ) , SLO ( I ) , STDEV (  I  ) 
WRITE(6,503) 

WRITE (6, 500)  (TLA( I), TLO (I), 1=1, NTAR) 

WRITE  (6,501)  ( SLA(  I)  ,SLO( I ) ,STDEV(  I)  ,1  =  1, NST) 
C 

DO  10  1=1, NST 
C 

C       COMPUTE    HEADING 
C 

IF( I.EQ.NST)    GOTO    5 

A=SLA( I+l)-SLA( I ) 

TH=(SLA(  I  +  D+SLACI  )  )/2. 

B=(SLO( 1+1 )-SLO( I ) )*COS(TH*RADDEG) 

IF( ABS( BJ.LT.l .E-6)     G0T03 

FEAD=PI/2.-ATAN( A/B) 

IF(B.LT.O)     HEAD=HEAD+PI 

GOTO    5 
3  HEAD=0. 

IFU.LT.O.)  HEAD  =  PI 
5    HD(I)=HEAD/RADDEG 
C 

C   COMPUTE  BEARING  ANGLES  AND  DIST  FROM  I-TH  STATION  TO  J-TH 
C 

DO  10  J=1,NTAR 

IFLAG=2 

IF(SLO( I).GT.TLO(J) )  IFLAG=1 

CA=ABS( SLO(  I  )-TLO( J) )*RADDEG 

SB=(90.-TLA( J) )*RADDEG 

SC=(90.-SLA( I ) )*RADDEG 

SA=ARCOS(COS(SB)*COS( SC ) +S IN( SB )*SIN< SC )*COS( CA) ) 

CB=ARCOS( (COS(SS)-COS(SA)*COS(SC) )/(SIN(SA)*SIN(SC))) 

THT( I , J)=CB/RADDEG 

GOTO  (4,7),IFLAG 
4   THT( I, J)=360.-THT( I ,J) 
7  CONTINUE 
C 

C   DIST  CALCULATION 
C 

DIST( I , J)=SA*RADIUS 
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C   STDEV  AS  A  FUNCTION  OF  REL.  BEARING 
C 

ANG=HEAD-THT(I ,J )*RADDEG 

STEEST( I ,J)=STDEV( I )* ( COS ( ANG )**2+l . ) 
10  CONTINUE 
C 

C   OUTPUT  RESULTS 
C 

DO  54  K  =  l,4 

IL=K*6 

Il=IL-5 

WRITE  (6,502) 

DO  53  1  =  11, IL 

IF(I.GT.NST)   GOTO  54 

WRITE  (6,253)  I 

53  WRITE  (6,255)  ( J ,THT ( I , J ) , DI ST (  I ,  J  )  , J  =  l , NTAR ) 

54  CONTINUE 
C 

C   COMPUTE  X-Y  POSITIONS  OF  STATIONS 
C 

DO  1000  1=1, NST 

CALL  SHIFT(RADIUS,SLO( I ) , SL AC  I ) , X ( I ) , Y ( I ) ) 

XX(I )=X( I) 

YY( I)=Y( I) 

1000  CONTINUE 

CALL  DRAW(NST,X,Y,1,0,LB, ITL,.7,.7,0,6,2,2,6,6,0,L) 

CALL  DRAW(NST,X,Y,2,1,LBL,ITL,.1,.05,0,9,2,2,9,10,0,L) 
C 

C   COMPUTE  X-Y  POSITIONS  OF  TARGETS 
C 

DO  1001  1=1, NTAR 

CALL  SHIFT(RADIUS,TLO( I ) ,TLA ( I ) , X ( I ) , Y( I  )  ) 

1001  CONTINUE 

CALL  DRAW(NTAR,X,Y,2,2,LBL,ITL,.1,.1,0,9,2,2,9,10,0,L) 
RETURN 
100  F0RMATU2) 
102  F0RMAT(3F10.4) 
253  F0RMAT(///14X, ' BEARING  AND  DISTANCE  FROM  STATION', 13, 

£•  TO  TARGET  INDICATED1,//) 
255  FORMAT( I 16,2F20.2) 

500  FORMAT (///14X, ' TARLAT' ,4X , • T ARLON • , // ( 1 2X, 2F 10 . 5 ) ) 
5  01  FORMAT (///14X,  • STAL AT ' , 4X , • ST ALON • , 5X , ■ STDEV  ,// 

£(12X,3F10.5)  ) 

502  FORMAT( 1H1 ) 

503  FORMAT( 1H1 ,////14X, 'COMPUTER  RESULTS  FOR  SIX  TARGETS' 
£,«  AND  50%  DETECTION'  ) 

END 

SUBROUTINE  CONVT  ( PHI ,TH1 ,0B 1 , A , B, C ) 
C 

DATA  PI/3.141593  /, RD/ . 0174533/ 

PH=PH1*RD 

TH=TH1*RD 

0B=0B1*R0 

ST=SIN(TH) 

CT=COS(TH) 

SP=SIN(PH) 

CP=COS( PH) 

CO=COS(OB) 

SO=SIN(OB) 

DX=-SP*C0*CT-ST*SO 

DY=SO*CT-SP*CO*ST 

0Z=CP*CO 

X1=CP*CT 

Yl=CP*ST 

A=Y1*DZ-SP*DY 

B=SP*DX-X1*DZ 

C=X1*DY-Y1*DX 

RETURN 

END 
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SUBROUTINE  CORREL  ( IC, NO, NT ,OBS , X , GATE, ADDGT , KK , NST A) 
C 

C  LIST  OF  VARIABLES 
C 

C  X(K?L)       TRACK  MATRIX 

C  OBS(M,N)     OBSERVATION  MATRIX 

C  C(I,J)       CORRELATION  MATRIX 
C 

C   NOTE  SUBSCRIPTS  ARE  USED  AS  INDICATED  ABOVE  IF  POSSIBLE 
C 

C  S(I)   DISTANCE  IN  ANGLE  OF  OBSERVATION  I  FROM  TRACK  J 

C  =101.  IF  UNCORRELATED 

C  S(J)   DISTANCE  IN  ANGLE  OF  OBSERVATION  I  FROM  TRACK  J 

C  =101.  IF  UNCORRELATED 
C 

c 

INTEGER  C(20,20) ,SUM 
COMMON  /A/  NTR,NJ,TH1,TD,NFRST 

DIMENSION  X(5,20),GATE(5) ,  OBS ( 5 , 20 ) , S ( 20 ) , ADDGT ( 20 ) , 
1KK(20) ,NJ(20,20) 
DIMENSION  TD( 2 0,20)  ,NFRST( 20) ,IC( 20) ,NTR(20) ,TH1(20) 
C 

DATA  TDM/2. 2918/, RD/. 0174533/ 
IF(NSTA.GT.2)  GOTO  11 
DO  10  1=1,20 
DO  10  J=l,20 

10  C(I,J)=0 

11  CONTINUE 
NMN=NT 
NMX=NT 

IF(NC.GT.NT)  NMX=NO 
IF(NO.LT.NT)  NMN=NO 
DO  12  1=1, NMX 

IC( I )=0 

12  C( I, I)=0 
GT=GATF(5) 

C 

WRITE  (6,2) 

WRITE  (6,9) 

WRITE  (6,3)  ( ( X(K,L) ,K=1,5) ,L=1,NT) 

WRITE  (6,5) 

WRITE  (6,9) 

WRITE (6, 3)  ( (OBS(M,N)  ,M=1,5)  ,N=1,N0) 

DO  18  L=1,NT 

GATE(5)=GT+ADDGT(L) 

DO  18  N=1,N0 

DO  17  K=l,5 

TEMP=OBS(K,N)-X(K,L) 

IF( ABS(TEMP) .GT.GATE(K) )  GOTO  18 
17   CONTINUE 

C(N,L)=1 
18  CONTINUE 

GATE(5)=GT 
C 

WRITE  (6,6) 

WRITE  (6,7)  (NTR( I) ,1=1,10) 

DO  180  1=1, NO 
180  WRITE  (6,4)  NJ ( NSTA , I ) , ( C( I , J ) , J= 1 ,NT ) 
C 

c 

C  RULE  1 
C 

C  A  TRACK  CORRELATING  WITH  SEVERAL  OBSERVATIONS  REJECTS 

C  ANY  OBSERVATION  HELD  IN  COMMON  WITH  AN  OTHER  TRACK,  IF 

C  THE  COMMON  OBSERVATION  IS  THE  ONLY  OBSERVATION 

C  CORRELATING  WITH  THE  OTHER  TRACK 
C 


DO  19  J=1,NT 
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c 

C      COMPUTE  #  OF  OBSERVATIONS  CORRELATED  TO  TRACK  J 
C 

20   SUM=0 

DO  25  1=1, NC 
25   SUM=C( I  ,J)+SUM 
C 

C      CHECK  IF  ONLY  1  OBSERVATION  CORRELATED  . 
C 

IF(SUM.LE.l)  GOTO  19 

DO  119  1=1, NO 

IF(C( I, J  ).E0.O )  GOTO  119 

DO  21  JP=1,NT 

IF(C( I, JP) .EQ.O)  GOTO  21 
C 

C      COMPUTE  #  OF  OBSERVATIONS  CORRELATED  TO  TRACK  JP 
C 

SUM=0 

DO  23  IP=1,N0 
23   SUM=SUM+C( IP, JP) 
C 

C      CHECK  IF  ONLY  1  OBSERVATION  CORRELATED 
C 

IF(SUM.NE.l)  GOTO  21 

C(I,J)=0 

GOTO  20 
21  CONTINUE 
119  CONTINUE 
19  CONTINUE 
C 
C 

C      RULE  2 
C 

C    A  TRACK  CORRELATING  WITH  SEVERAL  OBSERVATIONS,  SOME  OF 
C    WHICH  ARE"  NOT  HELD  IN  COMMON  WITH  OTHER  TRACKS,  REJECTS 
C    OBSERVATIONS  HELD  IN  CCMMON  WITH  OTHER  TRACKS 
C 

DO  31  J=1,NT 
C 

C      COMPUTE  #  OF  OBSERVATIONS  CORRELATED  TO  TRACK  J 
C 

SUM=0 

DO  32  1=1, NC 
32   SUM=SUM+C( I, J) 
C 

C      CHECK  IF  ONLY  1  OBSERVATION  CORRELATED 
C 

IF(SUM.LE.l)  GOTO  31 

DO  131  1=1, NO 

IF(C(  I, J). EQ.O)  GOTO  131 
C 

C      OBSERVATION  I  IS  ONE  OF  SEVERAL  OBSERVATIONS 
C      CORRELATED  TO  TRACK  J 

C      COMPUTE  #  OF  TRACKS  CORRELATED  TO  OBSERVATION  I 
C 

SUM=0 

DO  35  JP=1,NT 
35   SUM=SUM+C( I ,JP) 
C 

C      CHECK  WHETHER  IN  COMMON  WITH  OTHER  TRACKS 
C 

IF(SUM.GT.l)  GOTO  131 
C 

C      REJECT  OBSERVATIONS  IN  COMMON  WITH  OTHER  TRACKS 
C 

DO  231  K=1,N0 

IF(C(K, J). EQ.O)  GOTO  231 

SUM  =  0 

DO  39  JP=1,NT 
39   SUM=SUM+C(K,JP) 

IF(SUM.EQ.l)  GOTO  231 

C(K,J)=0 
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231  CONTINUE 

GOTO  31 
131  CONTINUE 
31  CONTINUE 
C 
C 

C      RULE  3 
C 

C    WHEN  SEVERAL  OBSERVATIONS  CORRELATE  WITH  ONE  TRACK,  THE 
C    CLOSEST  OBSERVATION  IS  CORRELATED  TO  THAT  TRACK 
C 

DO  41  J=1,NT 

SUM=0 
C 

C      COMPUTE  #    OF  TRACKS  CORRELATED. 
C 

DO  42  I=1,NC 
42   SUM=SUM+C( I,J) 

IF(SUM.LE.l )  GOTO  41 
C 

C      COMPUTE  DISTANCE 
C 

DO  44  1=1, NC 

S( I )=101.0 

IF(C( I, J) .EO.O)  GOTO  44 

S( I )=ABS(X( 5,J )-OBS(5,I ) ) 
44  CONTINUE 
C 

C      FIND  CLOSEST  OBSERVATION 
C 

NB=1 

DO  46  IB=2»NO 

IF(S(1) .LE.S( IB) )  GOTO  46 

S(1)=S( IB) 

NB=IB 
46  CONTINUE 

DO  49  1=1, NC 
49  C(  I, J)  =  0 
C 

C      RESET  PROPER  ELEMENT 
C 

C  (  NB  ,  J  )  =  1 
41  CONTINUE 
C 

c 

C      RULE  4 
C 

C    WHEN  SEVERAL  TRACKS  CORRELATE  WITH  ONE  OBSERVATION,  THE 

C    OBSERVATION  IS  CORRELATED  WITH  THE  CLOSEST  TRACK 

C 

c 
c 
c 


DO  51  1=1, NO 
C 

C      COMPUTE  #    OF  TRACKS  CORRELATED 


SUM=0 

DO  52  J=1,NT 
52   SUM=SUM+C( I, J) 

IF(SUM.LE.l )  GOTO  51 
C 

C      COMPUTE  DISTANCE 
C 

DO  54  J=1,NT 

S(J)=101.0 

IF(C(  I,  J)  .ECO)  GOTO  54 

S( J)=ABS(X(5,J )-OBS(5,I ) ) 
54  CONTINUE 
C 

C      FIND  CLOSEST  TRACK. 
C 

NC=1 

DO  56  JB=2,NT 

IF(S(1)  .LE.S(JB)  )  GOTO  56 
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S(1)=S(JB) 

NC  =  JB 
56  CONTINUE 

DO  59  J=1,NT 
59   C(I,J)=0 
C 
C      RESET  PROPER  ELEMENT 


C 


C(  I,NC)=1 
51  CONTINUE 


C 
C 

C      RULE  5 
C 

C    IF  AN  OBSERVATION  IS  FOUND  UNCORRELATED,  AND  AT  THE  SAME 
C    TIME  THERE  ARE  TRACKS  UNCORRELATED  TO  ANY  OBSERVATION, 
C    THEN  THIS  OBSERVATION  IS  CORRELATED  TO  THE  FIRST 
C    UNCORRELATED  TRACK  FOUND  WITHIN  THE  GATE 
C 

DO  515  1=1, NO 

DO  511  J=1,NT 

IF(C(  I, J  ).EQ.l )  GOTO  515 

511  CONTINUE 

DO  514  J=1,NT 
DO  512  K=l,NO 
IF(C(K, JJ.EG.l)  GOTO  514 

512  CONTINUE 
GATE(5)=GT+ADDGT( J ) 
DO  513  K=l,5 

TEMP  =OBS(K, I)-X(K, J) 

IF( ABS(TEMP) .GT.GATE(K) )  GOTO  514 

513  CONTINUE 
C(  I,J)  =  1 
GOTO  515 

514  CONTINUE 

515  CONTINUE 
GATE(5)=GT 

C 

WRITE  (6,8) 

WRITE  (6,7)  (NTRU  )  ,1  =  1,10) 

DO  510  1=1, NO 
510  WRITE  (6,4)  NJ ( NSTA , I ) ,  ( C ( I , J )  ,  J  =  l  ,NT  ) 
C 

C    REARRANGE  C-MATRIX  TO  DIAGONAL  FORM 
C 

L  =  0 
C 

C    SEARCH  DIAGONAL  FOR  ZEROES 
C 

DO  62  1=1, NC 

M=I 

IF(C(  I,  I  J.NE.l  )  GOTO  63 

IC(  I)=C(  I,  I  ) 

62  CONTINUE 
IF(NO.EQ.NT)  RETURN 
GOTO  75 

C 

C    SEARCH  M-TH  COLUMN  FOR  OFF  DIAGONAL  1 

C    INTERCHANGE  CORRESPONDING  ROWS 

C 

63  IF(M.EO.NO)  GOTO  66 
MM=M 

DO  65  M=MM,NO 

IF(M.GT.NMN)  GOTO  66 

MP=M+1 

DO  65  I=MP, NO 

IF(C( I,M).EQ.O)  GOTO  65 

C(  I,M)=0 

C(M,M)=1 

DO  64  J=MP,NT 

IF(M.GE.NT)  GOTO  640 
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IF (C( M,J).EQ.O )  GOTO  64 

C(MfJ)=0 

C(  I,J)=1 

GOTO  640 

64  CONTINUE 
640  CONTINUE 

TEMP=0BS(5,M) 
OBS(5,M)=OBS(5, I  ) 
0BS(5, I  )=TEMP 

65  CONTINUE 

66  CONTINUE 
C 

C    SEARCH  FOR  UNCORRELATED  OBSERVATIONS 


C 


DO  69  1=1, NC 

IF(C( I,  I  ).EQ.l)  GOTO   69 

IP=I+1 

IF(IP.GT.NT)  GOTO  675 

DO  67  J=IP,NT 

IF(C(  I, J  ).EQ.l )  GOTO   69 


67 

CONTINUE 

c 
c 
c 

675 

CONTINUE 

INITIALIZE    NEW    TRACK 

L=L-H 

NI=NT+L 

C(NI tNI )=1 

NFRST(NI )=NSTA 

KK(NI ) =NSTA-1 

NTR(NI  )=NJ(  N'STAt  I  ) 

THKNI  )=0BS(5,  I  ) 

DO    68    J  =  l,5 

X( J,NI )=OBS(J, I ) 

68 

CONTINUE 

KKNI=KK(NI ) 

TD(KKNI »NI )=TDM*SIN(X( 5,NI)*RD) 

69 

CONTINUE 

IF(L.EQ.O)     GOTO    700 

NTP=NT+1 

DO    70    I=NTP,NI 

DO    70    J  =  l,5 

70 

OBS( JTI  )=X( J , I) 

700 

CONTINUE 

c 

C  CONTINUE  TO  REARRANGE  THE  C-MATRIX 
C  SEARCH  I-TH  RCW  FOR  OFF  DIAGONAL  1 
C 

DO  74  1=1, NT 

IF(C(  I, I  ).EQ.l )  GOTO  74 

IP=I+1 

DO  73  J=IP, NT 

IF(C( I, JJ.EQ.O)  GOTO  73 

DO  71  K=J,NT 

IF(C( J, K).EQ.O)  GOTO  71 

C  (  I ,  K  )  =  1 

C( J,K)=0 

71  CONTINUE 
C(  I,  J)=0 
C  (  J  ,  J )  =  1 

DO  72  K=l,5 
TEMP=OBS(K, J) 
OBS(K, J)=OBS(K,I ) 

72  OBS(K, I )=TEMP 

73  CONTINUE 

74  CONTINUE 
C 

C    SEARCH  FOR  MISSING  DATA 

C    SET  KALMAN  FILTER  GAIN  INDEX 


C 


75   NT=NT+L 

DO  77  1=1, NT 
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IC( I)  =  C( 1,1 ) 

IF( IC( I ) .EQ.O)  GOTO  76 

X(4,I)=10 

GOTO  77 

76  KK(  I)=KK(I  )+l 
X(4,I )=X(4,I)-1. 

IF(KK( I ) .GE.NSTA)  KK(I)=NSTA-1 

77  CONTINUE 
RETURN 

2  FORMAT  (///14X, 'ORIGINAL  MATRIX  X  OF  PREDICTED  TRACKS' 

3  FORMAT  (/, (12X,5F10.3)  ) 

4  FORMAT  (121, 5X, 2013) 

5  FORMAT  (//, 14X, 18HINPUT  OBSERVATIONS) 

6  FORMAT  (//,14X, • INITIAL  CORRELATION  MATRIX'/) 

7  FORMAT  (/30X, 'TRACKS' , //26X, 101 3, /14X, 'OBSERVATIONS'/) 

8  FORMAT  (/14X,24HFINAL  CORRELATION  MATRIX,/) 

9  FORMAT  (/14X,'  FREQU' ,6X, ' PRF •  ,9X,  ' PW ,7X,  'CCUNT     DOA 
END 

SUBROUTINE  FIND  ( TLA,TLO , S IG, N, TH , SLA, SLO ) 

DIMENSION  A(20),B(20),C(20) ,XS(3) ,XM(3) , XI (200) ,X2( 200 
1X3(200) , SLA (20) , SLO (20) ,TH(20) 

DATA  RD/. 01744533/, R/3959./ 
NM=N-1 
M  =  0 

DO  1  J=1,N 
1    CALL  CONVT  ( SL A ( J ) , SLO( J ) ,TH ( J ) , A ( J ) , B( J ) , C ( J ) ) 
DO  3  0  1=1, NM 

k=hi 

DO  10  J=K,N 

M=M+1 

X3(M)=1. 

D=A( I)*B(J)-A( J)*B(  I) 

X1(M)  =  (B(I)*C(  J)-B(  J)*CU)  )/D 

X2CM)=( A(J)*C( I)-A(I)*C( J) )/D 

E=R/S0RT(X1(M) **2+X2(M)**2+l. ) 

X1(M)=X1(M)*E 

X2(M)=X2(M)*E 

10  X3(M)=X3(M)*E 
DO  11  1=1,3 
XM(I )=0. 

11  XS(I)=0. 

DO  12  1=1, M 
XM(1)=XM(1)  +  X1  (  I) 
XM(2)=XM(2)  +  X2 (  I  ) 

12  XM(3)=XM(3)+X3 (  I  ) 
IF(N.EQ.2)  GOTO  16 
DO  13  J=l,3 

13  XM( J)=XM( J )/M 
DO  14  1=1, M 

XS(1)=XS(1)+(X1(  I)-XM(l)  )**2 
XS(2)=XS(2)+(X2(I)-XM(2) )**2 

14  XS(3)=XS(3)+(X3( I )-XM( 3) )**2 
DO  15  J=l,3 

15  XS( J)=XS(J)/M 
SIG=SQRT(XS(1)+XS(2)+XS(3)  ) 

16  CONTINUE 
TL0=ATAN2(XM(2) ,XM( 1) ) /RD 
SQ=SQRT(XM( 1)**2+XM<2)**2) 
TLA=ATAN(XM(3)/SQ)/RD 
RETURN 

END 


78 


BIBLIOGRAPHY 


[1]   Kalman,  R.  E.,  "A  New  Approach  to  Linear  Filtering 

and  Prediction  Problems,"  Transactions  ASME  Journal 
Basic  Engineering,  Ser.  D,  vol.  82,  no.  2,  p.  35-45, 
March  19  60. 

[2]   Rauch,  H.  E. ,  "Solutions  to  the  Linear  Smoothing 
Problem, "  IEEE  Transactions  on  Automatic  Control, 
vol.  AC-8,  no.  4,  p.  371-372,  October  1963. 

[3]   Kanyusk,  A.  J.,  Singer,  R.  A.,  "Correlation  of  Multiple- 
Site  Track  Data,"  IEEE  Transaction  on  Aerospace  and 
Electronic  Systems,  vol.  AES-6,  no.  2,  p.  180-187, 
March  19  70. 

[4]  Titus,  H.  A.,  Pope,  W.  R. ,  "Multiple  Emitter  Airborne 
Direction-Finding  with  EOB  Utilization, "  Proc.  of  the 
Sixteenth  Annual  Joint  Electronic  Warfare  Conference, 
Naval  Postgraduate  School,  Monterey,  CAlif.,  Dec.  19  70. 


79 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  9  3940 

3.  Associate  Professor  H.  A.  Titus,  Code  52Ts        20 
Electrical  Engineering  Department 

Naval  Postgraduate  School 
Monterey,  California  9  3940 

4.  Associate  Professor  D.  E.  Kirk,  Code  52Ki         1 
Electrical  Engineering  Department 

Naval  Postgraduate  School 
Monterey,  California  9  3940 

5.  Felix  Pfendtner  1 
752  Bruchsal 

Michelsbergstr.  40,  GERMANY 

6.  Inspektion  des  Erziehungs-und  1 
Bildungsweseus  der  Marine 

294  Wilhelmshaven,  GERMANY 

7.  Dokumentationszentrale  der  Bundeswehr  [See]        i 
52  Bonn 

Friedrich  Ebert  Allee  34,  GERMANY 

8.  CDR  D.  Prentiss,  USN  1 
PM-7-21 

Naval  Material  Command 
Washington,  D.C.   20360 

9.  LCDR  R.  E.  Hudson,  USN  1 
AIR-53364C 

Naval  Air  Systems  Command 
Washington,  D.  C.   20  360 

10.  Mr.  James  W.  Hannan  1 
Airtronics,  Inc. 

SWL  Division 

Dulles  International  Airport 

Washington,  D.C.   20041 


80 


Distribution  List  (cont.)  No.  Copies 

11.  Commanding  Officer  1 
Naval  Missile  Center 

Naval  Air  Station 

Pt.  Mugu,  California   93041 

Attn:   Code  5341  (Van  Cleave) 

12.  Commanding  Officer  1 
Naval  Avionics  Facility 

21st  and  Arlington  Avenues 
Indianapolis,  Indiana   46218 
Attn:   Code  A20  (Brienza) 

13.  Mr.  Robert  Welsh  1 
Sperry  Rand  Corporation 

UNIVAC  Federal  Systems  Division 
UNIVAC  Park,  P.O.  Box  3525 
St.  Paul,  Minnesota   55101 

14.  Mr.  Lloyd  Trace  1 
United  Technology  Laboratories 

P.O.  Box  30102 
Dallas,  Texas   75230 

15.  Mr.  L.  J.  Polansky  1 
Lockheed  California  Company 

P.O.  Box  551 

Burbank,  California   91503 

16.  Dean  of  Research  Administration  1 
Code  02  3 

Naval  Postgraduate  School 
Monterey,  California   9  3940 


81 


82 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 

(Security  classification  ol  title,   body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  Is  classitied) 


I      originating    activity   (Corporate  author) 

Naval  Postgraduate  School 
Monterey,  California  9  3940 


2a.   REPORT    SECURITY    CLASSIFICATION 


Unclassified 


2b.     GROUP 


3      REPORT     TITLE 


Adaptive  Angle  Tracking  and  Correlation  for  Airborne  Direction-Finding 


4     DESCRIPTIVE   NOTES  (Type  of  report  and. inclusive  dates) 


Electrical  Engineer's  Thesis;   March  19  71 


5     au  THORiSi  (First  name,  middle  initial,  last  name) 


Felix   Pfendtner 


6      REPOR  T    DATE 

March    19  71 


7a.     TOTAL    NO.    OF    PAGES 

84 


7b.    NO.    OF    RE  FS 

4 


ta      CONTRACT    OR    GRANT    NO. 


6.    PROJEC  T    NO 


9o.    ORIGINATOR"S    REPORT    NUMBER(S) 


9b.    OTHEf*   REPORT   NOISI  (Any  other  numbers   that  may  be  assigned 
this  report) 


10      DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited 


II      SUPPLEMENTARY    NOTES 


12.    SPONSORING    MILITARY    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California  93940 


13.  ABSTRACT 


An  airborne  direction-finding  technique  capable  of  handling 
multiple  emitters  was  developed.   An  adaptive  gate  size  was  introduced 
in  the  track  correlations  based  on  covariance  relations  of  established 
tracks  and  observations.   Track  files  were  generated  based  upon  the 
existence  of  various  target  parameters;  i.e.,  frequency,  pulse  repeti- 
tion frequency,  pulse  width  and  direction-of-arrival .   To  test  the 
angular  resolution  capability  of  the  filter,  emitters  in  close  prox- 
imity to  each  other  with  identical  electronic  characteristics  were  used 
in  the  simulation.   Target  locations  are  calculated  in  a  cartesian 
coordinate  system  where  the  sphericity  of  the  earth  is  taken  into 
account,  and  with  appropriate  coordinate  transformation  computational 
simplicity  is  preserved. 


DD.f, 


nov  es  I  *v  /  O 

S/N    0101 -807-681 1 


(PAGE     1  ) 


UNCLASSIFIED 


Security  Classification 


A-31408 


UNCLASSIFIED 


Security  Classification 


KEY     WO  R  OS 


Angle  Tracking 
Kalman  Filter 
Direction-Finding 
Track  Correlation 


FORM 

I    NOV    6S 

S/N     0101-807-682  1 


1473     (BACK) 


84 


UNCLASSIFIED 


Security  Classification 


A-  3 1 409 


Th 


126277 


Thesi  s 
P4617 

c.l 


pfendtner 

Adaptive  angle 
tracking  and  cor- 
relation for  air- 
borne direction- 
finding- 


in 


4  ftS72 

OCT  31  65 

3      APR    86 

23  SfP   ^ 


o  0  9  6  I' 
2  5  9  7  tt 
5  0368 

302  17 
3697^ 


91  e 

cor- 

fir- 
n~ 


5fl 


f 


12527 


Pfendtner 

Adaptive   angle 
tracki  ng  and  cor- 
relation   for  ai  r- 
borne  di  rection- 
f i  ndi  ng. 


,hesP46i7  .rarkina  and  correlation 

Adaptive  angle  tr a9.K.' "?..„.,.«, u.» mil uiu Ml HH 


3  2768  001001910 

DUDLEY  KNOX  LIBRARY 


